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We present a measurement of the W boson mass using data collected by the D0 experiment at 
the Fermilab Tevatron during 1994-1995. We identify W bosons by their decays to ev final states. 
We extract the W mass, Mw, by fitting the transverse mass and transverse electron momentum 
spectra from a sample of 28,323 W ^ eu decay candidates. We use a sample of 3,563 dielectron 
events, mostly due to Z — > ee decays, to constrain our model of the detector response. From the 
transverse mass fit we measure Mw = 80.44 ± Q.lQ[stat) ± Q.Q7{syst) GeV. Combining this with 
our previously published result from data taken in 1992-1993, we obtain Mw = 80.43 ± 0.11 GeV. 
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I. INTRODUCTION 

In this article we describe the most precise measure- 
ment to date of the mass of the W boson, using data 
collected in 1994-1995 with the D0 detector at the Fer- 
milab Tevatron pp collider ■ 

The study of the properties of the W boson began 
in 1983 with its discovery by the UAl |§ and UA2 |] 
collaborations at the CERN pp collider. Together with 
the discovery of the Z boson in the same year it 
provided a direct confirmation of the unified model of 
the weak and electromagnetic interactions which - 
together with QCD - is now called the Standard Model. 

Since the W and Z bosons are carriers of the weak 
force, their properties are intimately coupled to the struc- 
ture of the model. The properties of the Z boson have 
been studied in great detail in e+e~ collisions The 
study of the W boson has proven to be significantly more 
difficult, since it is charged and therefore can not be res- 
onantly produced in e~^e~ collisions. Until recently its 
direct study has therefore been the realm of experiments 
at pp colliders which have performed the most precise 
direct measurements of the W boson mass ]l0|-|l2[. Di- 
rect measurements of the W boson mass have also been 
carried out at LEP2 p"^-p^ using nonresonant W pair 
production. A summary of these measurements can be 



found in Table XV at the end of this article. 



The Standard Model links the W boson mass to other 
parameters. 



Mw 
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In the "on shell" scheme 1^11 
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where 9\y is the weak mixing angle. Aside from the radia- 
tive corrections Ar, the W boson mass is thus determined 
by three precisely measured quantities, the mass of the Z 
boson Mz 1^, the Fermi constant Gp (iS) and the elec- 
tromagnetic coupling constant a evaluated at = 



Mz = 91.1865 ± 0.0020 GeV, (3) 
Gf = (1.16639 ± 0.00002) x 10"'^ GeV"^ (4) 
a = (128.896 ±0.090)-^ (5) 

From the measured W boson mass we can derive the size 
of the radiative corrections Ar. Within the framework 
of the Standard Model, these corrections are dominated 
by loops involving the top quark and the Higgs boson 
(see Fig. The correction from the tb loop is substan- 
tial because of the large mass difference between the two 
quarks. It is proportional to for large values of the 
top quark mass mj. Since mj has been measured [ pO[ , 
this contribution can be calculated within the Standard 



Model. For a large Higgs boson mass, mn, the correc- 
tion from the Higgs loop is proportional to InniH- In ex- 
tensions to the Standard Model new particles may give 
rise to additional corrections to the value of Mw- In 
the Minimal Supersymmetric extension of the Standard 
Model (MSSM), for example, additional corrections can 
increase the predicted W mass by up to 250 MeV pT[ |. 

A measurement of the W boson mass therefore con- 
stitutes a test of the Standard Model. In conjunction 
with a measurement of the top quark mass the Standard 
Model predicts Mw up to a 200 MeV uncertainty due to 
the unknown Higgs boson mass. By comparing with the 
measured value of the W boson mass we can constrain 
the mass of the Higgs boson, the agent of the electroweak 
symmetry breaking that has up to now eluded experi- 
mental detection. A discrepancy with the range allowed 
by the Standard Model could indicate new physics. The 
experimental challenge is thus to measure the W boson 
mass to sufficient precision, about 0.1%, to be sensitive 
to these corrections. 



II. OVERVIEW 

A. Conventions 

We use a Cartesian coordinate system with the z-axis 
defined by the direction of the proton beam, the x-axis 
pointing radially out of the Tevatron ring and the y-axis 
pointing up. A vector p is then defined in terms of its 
projections on these three axes, Px, Py, Pz- Since pro- 
tons and antiprotons in the Tevatron are unpolarized, 
all physical processes are invariant with respect to ro- 
tations around the beam direction. It is therefore con- 
venient to use a cylindrical coordinate system, in which 
the same vector is given by the magnitude of its compo- 
nent transverse to the beam direction, px, its azimuth 0, 
and Pz- In pp collisions the center of mass frame of the 
parton-parton collisions is approximately at rest in the 
plane transverse to the beam direction but has an unde- 
termined motion along the beam direction. Therefore the 
plane transverse to the beam direction is of special im- 
portance and sometimes we work with two-dimensional 
vectors defined in the x-y plane. They are written with a 
subscript T, e.g. px- We also use spherical coordinates by 
replacing pz with the colatitude 9 or the pseudorapidity 
rj = — In tan (0/2). The origin of the coordinate system 
is in general the reconstructed position of the pp interac- 
tion when describing the interaction, and the geometrical 
center of the detector when describing the detector. For 
convenience, we use units in which c = fi = 1. 



B. W and Z Boson Production and Decay 

In pp collisions at y/s = 1.8 TeV, W and Z bosons are 
produced predominantly through quark- antiquark anni- 
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hilation. Figure || shows the lowest-order diagrams. The 
quarks in the initial state may radiate gluons which are 
usually very soft but may sometimes be energetic enough 
to give rise to hadron jets in the detector. In the reac- 
tion the initial proton and antiproton break up and the 
fragments hadronize. We refer to everything except the 
vector boson and its decay products collectively as the 
underlying event. Since the initial proton and antiproton 
momentum vectors add to zero, the same must be true 
for the vector sum of all final state momenta and there- 
fore the vector boson recoils against all particles in the 
underlying event. The sum of the transverse momenta of 
the recoiling particles must balance the transverse mo- 
mentum of the boson, which is typically small compared 
to its mass but has a long tail to large values. 

We identify W and Z bosons by their leptonic decays. 
The D0 detector (Sec. [II) is best suited for a precision 
measurement of electrons and positronsj^, and we there- 
fore use the decay channel W ^ eiy to measure the W 
boson mass. Z — > ee decays serve as an important cal- 
ibration sample. About 11% of the W bosons decay to 
ev and about 3.3% of the Z bosons decay to ee. The 
leptons typically have transverse momenta of about half 
the mass of the decaying boson and are well isolated from 
other large energy deposits in the calorimeter. Interme- 
diate vector boson decays are the dominant source of 
isolated high-p^ leptons at the Tevatron, and therefore 
these decays allow us to select a clean sample of W and 
Z boson decays. 



C. Event Characteristics 

In events due to the process pp — > (W ev) + X, 
where X stands for the underlying event, we detect the 
electron and all particles recoiling against the W with 
pseudorapidity —4 < 77 < 4. The neutrino escapes unde- 
tected. In the calorimeter we cannot resolve individual 
recoil particles, but we measure their energies summed 
over detector segments. Recoil particles with |?7| > 4 
escape unmeasured through the beampipe, possibly car- 
rying away substantial momentum along the beam direc- 
tion. This means that we cannot measure the sum of the 
^-components of the recoil momenta, Uz , precisely. Since 
these particles escape at a very small angle with respect 
to the beam, their transverse momenta are typically small 
and can be neglected in the sum of the transverse recoil 
momenta, ut- We measure ut by summing the observed 
energy flow vectorially over all detector segments. Thus, 
we reduce the reconstruction of every candidate event to 
a measurement of the electron momentum p[e) and ut- 



Since the neutrino escapes undetected, the sum of all 
measured final state transverse momenta does not add to 
zero. The missing transverse momentum j^t, required to 
balance the transverse momentum sum, is a measure of 
the transverse momentum of the neutrino. The neutrino 
momentum component along the beam direction cannot 
be determined, because Uz is not measured well. The 
signature ofaW ev decay is therefore an isolated high- 
Pt electron and large missing transverse momentum. 

In the case oi Z ^ ee decays the signature consists 
of two isolated high-pT electrons and we measure the 
momenta of both leptons, p{ei) and p{e2), and ut in 
the detector. 



D. Mass Measurement Strategy 

Since Pz{v) is unknown, we cannot reconstruct the ep 
invariant mass for W —f ev candidate events and there- 
fore must resort to other kinematic variables for the mass 
measurement. 

For recent measurements [[U)H12| the transverse mass. 



rriT 



^2pT{e)pT{v) (1 - cos (</)(e) - ^{v))) (6) 



was used. This variable has the advantage that its spec- 
trum is relatively insensitive to the production dynamics 
of the W . Corrections to rriT due to the motion of the W 
are of order (qr/Mn/)^, where qt is the transverse mo- 
mentum of the W boson. It is also insensitive to selection 



^In the following we use "electron" generically for both elec- 
trons and positrons. 



biases that prefer certain event topologies (Sec. VIC). 
However, it makes use of the inferred neutrino pt and is 
therefore sensitive to the response of the detector to the 
recoil particles. 

The electron pT spectrum provides an alternative mea- 
surement of the W mass. It is measured with better reso- 
lution than the neutrino pT and is insensitive to the recoil 
momentum measurement. However, its shape is sensitive 
to the motion of the W and receives corrections of order 
qx/Mw. It thus requires a better understanding of the 
W boson production dynamics than the ttit spectrum. 

The niT and prie) spectra thus provide us with two 
complementary measurements. This is illustrated in 
Figs. ^ and |^, which show the effect of the motion of the 
W bosons and the detector resolutions on the shape of 
each of the two spectra. The solid line shows the shape of 
the distribution before the detector simulation and with 
qT=0- The points show the shape after qr is added to 
the system, and the shaded histogram also includes the 
detector simulation. We observe that the shape of the 
TTiT spectrum is dominated by detector resolutions and 
the shape of the pt{s) spectrum by the motion of the 
W. By performing the measurement using both spectra 
we provide a powerful cross-check with complementary 
systematics. 

Both spectra are equally sensitive to the electron en- 
ergy response of the detector. We calibrate this re- 
sponse by forcing the observed dielectron mass peak in 
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the Z ^ ee sample to agree with the known Z mass ^ 
(Sec. V]). This means that we effectively measure the 



ratio of W and Z masses, which is equivalent to a mea- 
surement of the W mass because the Z mass is known 
precisely. 

To carry out these measurements we perform a maxi- 
mum likelihood fit to the spectra. Since the shape of the 
spectra, including all the experimental effects, cannot be 
computed analytically, we need a Monte Carlo simulation 
program that can predict the shape of the spectra as a 
function of the W mass. To perform a measurement of 
the W mass to a precision of order 100 MeV we have to 
estimate individual systematic effects to 10 MeV. This 
requires a Monte Carlo sample of 2.5 million accepted W 
bosons for each such effect. The program therefore must 
be capable of generating large samples in a reasonable 
time. We achieve the required performance by employ- 
ing a parameterized model of the detector response. 

We next summarize the aspects of the accelerator 
and detector that are important for our measurement 
(Sec. ni). Then we describe the data selection (Sec. IV) 
and the fast Monte Carlo model (Sec. Most param- 
eters in the model are determined from our data. We 
describe the determination of the various components of 
the Monte Carlo model in Sees. [v|-px[ After tuning 
the model we fit the kinematic spectra (Sec. per- 
form some consistency checks (Sec . pC^ , and discuss the 
systematic uncertainties (Sec. KII). Section [XIIl sum- 
marizes the results and presents the conclusions. 



III. EXPERIMENTAL SETUP 
A. Accelerator 

The Fermilab Tevatron ||2^ collides proton and an- 
tiproton beams at a center-of-mass energy of -^/s = 1.8 
TeV. Six bunches each of protons and antiprotons cir- 
culate around the ring in opposite directions. Bunches 
cross at the intersection regions every 3.5 ^s. During 
the 1994-1995 running period, the accelerator reached a 
peak luminosity of 2.5 x 10^^cm~^s~^ and delivered an 
integrated luminosity of about 100 pb~^. 

The Tevatron tunnel also houses a 150 GeV proton 
synchrotron, called the Main Ring, which is used as an 
injector for the Tevatron. The Main Ring also serves 
to accelerate protons for antiproton production during 
collider operation. Since the Main Ring beampipe passes 
through the outer section of the D0 calorimeter, passing 
proton bunches give rise to backgrounds in the detector. 
We eliminate this background using timing cuts based on 
the accelerator clock signal. 



B. Detector 



1. Overview 

The D0 detector consists of three major subsystems: 
a central detector, a calorimeter (Fig. and a muon 
spectrometer. It is described in detail in Ref. |23|. We 
describe only the features that are most important for 
this measurement. 



2. Central Detector 

The central detector is designed to measure the trajec- 
tories of charged particles. It consists of a vertex drift 
chamber, a transition radiation detector, a central drift 
chamber (CDC), and two forward drift chambers (FDC). 
There is no central magnetic field. The CDC covers the 
region |?7| < 1.0. It is a jet-type drift chamber with delay 
lines to give the hit coordinates in the r-z plane. The 
FDC covers the region 1.4 < I77I < 3.0. 



3. Calorimeter 

The calorimeter is the most important part of the de- 
tector for this measurement. It is a sampling calorime- 
ter and uses uranium absorber plates and liquid argon 
as the active medium. It is divided into three parts: a 
central calorimeter (CC) and two end calorimeters (EC), 
each housed in its own cryostat. Each is segmented into 
an electromagnetic (EM) section, a fine hadronic (FH) 
section, and a coarse hadronic (CH) section, with in- 
creasingly coarser sampling. The CC-EM section is con- 
structed of 32 azimuthal modules. The entire calorimeter 
is divided into about 5000 pseudo-projective towers, each 
covering 0.1x0.1 in 77 x The EM section is segmented 
into four layers, 2, 2, 7, and 10 radiation lengths thick. 
The third layer, in which electromagnetic showers typ- 
ically reach their maximum, is transversely segmented 
into cells covering 0.05x0.05 in 77 x 0. The hadronic 
section is segmented into four layers (CC) or five lay- 
ers (EC). The entire calorimeter is 7-9 nuclear interac- 
tion lengths thick. There are no projective cracks in the 
calorimeter and it provides hermetic and almost uniform 
coverage for particles with \r]\ < 4. Figure]^ shows a view 
of the calorimeter and the central detector. 

The signals from arrays of 2 x 2 calorimeter towers, cov- 
ering 0.2x0.2 in 77 X (/), are added together electronically 
for the EM section only and for all sections, and shaped 
with a fast rise time for use in the Level 1 trigger. We 
refer to these arrays of 2 x 2 calorimeter towers as "trigger 
towers" . 

Figure ^ shows the pedestal spectrum of a calorimeter 
cell. The spectrum has an asymmetric tail from ioniza- 
tion caused by the intrinsic radioactivity of the uranium 
absorber plates. The data are corrected such that the 
mean pedestal is zero for each cell. To reduce the amount 
of data that have to be stored, the calorimeter readout 
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is zero-suppressed. Only cells with a signal that deviates 
from zero by more than twice the rms of the pedestal 
distribution are read out. This region of the pedestal 
spectrum is indicated by the shaded region in Fig. |^. 
Due to its asymmetry, the spectrum does not average to 
zero after zero-suppression. Thus the zero-suppression 
effectively causes a pedestal shift. 

The liquid argon has unit gain and therefore the 
calorimeter response was extremely stable during the en- 
tire run. Figure |^ shows the response of the liquid argon 
as monitored with radioactive sources of a and (3 parti- 
cles. Figures || and ^ show the gains and pedestals of a 
typical readout channel throughout the run. 

The EM calorimeter provides a measurement of energy 
and position of the electrons from the W and Z decays. 
Due to the fine segmentation of the third layer, we can 
measure the position of the shower centroid with a pre- 
cision of 2.5 mm in the azimuthal direction and 1 cm in 
the z-direction. 

We study the response of the EM calorimeter to elec- 
trons in beam tests |Q . To reconstruct the electron en- 
ergy we add the signals Ui observed in each EM layer 
{i = 1 ... 4) and the first FH layer (i = 5) of an array of 
5x5 calorimeter towers, centered on the most energetic 
tower, weighted by a layer dependent sampling weight Si, 



E : 



To determine the sampling weights we minimize 



2 (Pbcam 



E) 



a 



(7) 



(8) 



EM 



where the sum runs over all events and item is the reso- 
lution given in Eq. |. We obtain A ^ 2.96 MeV/ADC 
count, Sem = -347 MeV, si = 1.31, S2 = 0.85, S4 = 0.98, 
and S5 = 1.84. We arbitrarily fix S3 = 1. The value of 
Sem depends on the amount of dead material in front of 
the calorimeter. The parameters si to S4 weight the four 
EM layers and 55 the first FH layer. Figure 10 shows 



the fractional deviation of as a function of the beam 
momentum Pboam- Above 10 GeV they deviate by less 
than 0.3% from each other. 

The fractional energy resolution can be parameterized 
as a function of electron energy using constant, sampling, 
and noise terms as 



^EM 



(9) 



with CEM = 0.003, SEM = 0.135 GeVi/2 |2|,|2^, and 
"■EM = 0.43 GeV in the central calorimeter. The angle 9 
is the colatitude of the electron. Figure |ll| shows the frac- 
tional electron energy resolution versus beam momentum 
for a CC-EM module. The line shows the parametriza- 
tion of the resolution from Eq. ^ 



4. Luminosity Monitor 

Two arrays of scintillator hodoscopes, mounted in front 
of the EC cryostats, register hits with a 220 ps time reso- 
lution. They serve to detect that an inelastic pp interac- 
tion has taken place. The particles from the breakup of 
the proton give rise to hits in the hodoscopes on one side 
of the detector that are tightly clustered in time. The de- 
tector has a 91% acceptance for inelastic pp interactions. 
For events with a single interaction the location of the 
interaction vertex can be determined with a resolution of 
3 cm from the time difference between the hits on the two 
sides of the detector for use in the Level 2 trigger. This 
array is also called the Level trigger because the detec- 
tion of an inelastic pp interaction is a basic requirement 
of most trigger conditions. 



5. Trigger 

Readout of the detector is controlled by a two-level 
trigger system. 

Level 1 consists of an and-or network, that can be pro- 
grammed to trigger on a pp crossing if a number of prese- 
lected conditions are true. The Level 1 trigger decision is 
taken within the 3.5 /is time interval between crossings. 
As an extension to Level 1, a trigger processor (Level 1.5) 
may be invoked to execute simple algorithms on the lim- 
ited information available at the time of a Level 1 accept. 
For electrons, the processor uses the energy deposits in 
each trigger tower as inputs. The detector cannot ac- 
cept any triggers until the Level 1.5 processor completes 
execution and accepts or rejects the event. 

Level 2 of the trigger consists of a farm of 48 VAXsta- 
tion 4000's. At this level the complete event is available. 
More sophisticated algorithms refine the trigger decisions 
and events are accepted based on preprogrammed condi- 
tions. Events accepted by Level 2 are written to magnetic 
tape for offline reconstruction. 



IV. DATA SELECTION 
A. Trigger 

The conditions required at trigger Level 1 for W and 
Z candidates are: 

• pp interaction: Level hodoscopes register hits con- 
sistent with a pp interaction. This condition ac- 
cepts 98.6% of all W and Z bosons produced. 

• Main Ring Veto: No Main Ring proton bunch 
passes through the detector less than 800 ns be- 
fore or after the pp crossing and no protons were 
injected into the Main Ring less than 400 ms be- 
fore the pp crossing. 
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• EM trigger towers: There are one or more EM trig- 
ger towers with E sin 6 > T, where E is the energy 
measured in the tower, 6 its angle with the beam 
measured from the center of the detector, and T a 
programmable threshold. This requirement is fully 
efficient for electrons with px > 2T. 

The Level 1.5 processor recomputes the transverse elec- 
tron energy by adding the adjacent EM trigger tower with 
the largest signal to the EM trigger tower that exceeded 
the Level 1 threshold. In addition, the signal in the EM 
trigger tower that exceeded the Level 1 threshold must 
constitute at least 85% of the signal registered in this 
tower if the hadronic layers are also included. This EM 
fraction requirement is fully efficient for electron candi- 
dates that pass our offline selection (Sec. IV D| ). 

Level 2 uses the EM trigger tower that exceeded the 
Level 1 threshold as a starting point. The Level 2 al- 
gorithm finds the most energetic of the four calorimeter 
towers that make up the trigger tower, and sums the 
energy in the EM sections of a 3x3 array of calorime- 
ter towers around it. It checks the longitudinal shower 
shape by applying cuts on the fraction of the energy in 
the different EM layers. The transverse shower shape 
is characterized by the energy deposition pattern in the 
third EM layer. The difference between the energies in 
concentric regions covering 0.25x0.25 and 0.15x0.15 in 
T] X (j) must be consistent with an electron. Level 2 also 
imposes an isolation condition requiring 



Ei Ei sir 



PT 



< 0.15 



PT 



(10) 



where t he sum runs over all cells within a cone of radius 
R = \/ -I- Ary^ = 0.4 around the electron direction 
and Pt is the transverse momentum of the electron . 

The Pt of the electron computed at Level 2 is based 
on its energy and the z-position of the interaction ver- 
tex measured by the Level hodoscopes. Level 2 ac- 
cepts events that have a minimum number of EM clusters 
that satisfy the shape cuts and have pt above a prepro- 
grammed threshold. Figure ^ shows the measured rela- 
tive efficiency of the Level 2 electron filter versus electron 
Pt for a Level 2 pT threshold of 20 GcV. We determine 
this efficiency using Z data taken with a lower threshold 
value (16 GeV). The efficiency is the fraction of electrons 
above a Level 2 pT threshold of 20 GeV. The curve is the 
parameterization used in the fast Monte Carlo. 

Level 2 also computes the missing transverse momen- 
tum based on the energy registered in each calorimeter 
cell and the vertex z-position. We determine the effi- 
ciency curve for a 15 GeV Level 2 requirement from 
data taken without the Level 2 condition. Figure 
shows the measured efficiency versus pt{v)- The curve is 
the parameterization used in the fast Monte Carlo. 



B. Reconstruction 



1. Electron 

We identify electrons as clusters of adjacent calorime- 
ter cells with significant energy deposits. Only clusters 
with at least 90% of their energy in the EM section and at 
least 60% of their energy in the most energetic calorime- 
ter tower are considered as electron candidates. For most 
electrons we also reconstruct a track in the CDC or FDC 
that points towards the centroid of the cluster. 

We compute the electron energy E[e) from the signals 
in all cells of the EM layers and the first FH layer in 
a window covering 0.5x0.5 in 77 x and centered on the 
tower which registered the highest fraction of the electron 
energy. In the computation we use the sampling weights 
and calibration constants determined using the testbeam 
data (Sec. Ill B 3D except for the offs et (5e m; which we 
take from an in situ calibration (Sec. VI D), i.e. (5em = 
—0.16 GeV for electrons in the CC. 

The calorimeter shower centroid position (xcai, J/cai, 
Zeal), the center of gravity of the track (xtrk, 2/trk, ^trk) 
and the proton beam trajectory define the electron di- 
rection. The shower centroid algorithm is documented in 
Appendix The center of gravity of the CDC track 
is defined by the mean hit coordinates of all the de- 
lay line hits on the track. The calibration of the mea- 
sured z-coordinates contributes a significant systematic 
uncertainty to the W boson mass measurement and is 
described in Appendices ^ and Using tracks from 
many events reconstructed in the vertex drift chamber, 
we measure the beam trajectory for every run. The clos- 
est approach to the beam trajectory of the line through 
shower centroid and track center of gravity defines the 
position of the interaction vertex (xvtx, 2/vtx, -^vtx)- In 
Z ^ ee events we may have two electron candidates with 
tracks. In this case we take the point midway between 
the vertex positions determined from each electron as 
the interaction vertex. Using only the electron track to 
determine the position of the interaction vertex, rather 
than all tracks in the event, makes the resolution of this 
measurement less sensitive to the luminosity and avoids 
confusion between vertices in events with more than one 
pp interaction. 

We then define the azimuth 0(e) and the colatitude 
Q(e) of the electron using the vertex and the shower cen- 
troid positions. 



tan (/)(e) 



^/cal 2/vtx 
^cal ^vtx 



(11) 



tan0(e) = V^^'ai + ^cai V4tx + 2/vtx ^ ^^3) 



^cal ^vtx 



Neglecting the electron mass, the momentum of the elec- 
tron is given by 



sin6'(e) cos 0(e) 
p(e)=E(e)\ sin6'(e) sin(/)(e) 
cos 0(e) 



(13) 
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2. Recoil 

We reconstruct the transverse momentum of all par- 
ticles recoiling against the W oi Z boson by taking the 
vector sum 



ut 



Ei sin 6i 



cost 
sin f 



(14) 



where the sum runs over all calorimeter cells that were 
read out, except those that belong to electron clusters. 
Ei are the cell energies, and (pi and 9i are the azimuth 
and colatitude of the center of cell i with respect to the 
interaction vertex. 



3. Derived Quantities 



In the case of Z 
momentum 



p(ee) = p{ei) + p{e2 
and the dielectron invariant mass 



ee decays we define the dielectron 

(15) 



m{ee) = y/2E{ei)E{e2)il — cosw), 



(16) 



where uj is the opening angle between the two electrons. 
It is useful to define a coordinate system in the plane 
transverse to the beam that depends only on the electron 
directions. We follow the conventions first introduced by 
UA2 [|l0| and call the axis along the inner bisector of the 
two electrons the 77-axis and the axis perpendicular to 
that the ^-axis. Projections on these axes are denoted 
with subscripts 77 or ^. Figure |lj illustrates these defini- 
tions. 

In case oi W ^ ev decays we define the transverse 
neutrino momentum 



(17) 



and the transverse mass (Eq. Useful quantities are 
the projection of the transverse recoil momentum on the 
electron direction, 



'II 



UT -prie), 



(18) 



and the projection on the direction perpendicular to the 
electron direction. 



U_L = UT ■ {pT{e) X z). 



(19) 



Figure 15 illustrates these definitions. 



the calorimeter module boundaries in the CC-EM and 
from the edges of the calorimeter by cutting on Ac/) and 
Zeal- We also remove electrons for which the z-position of 
the track center of gravity is near the edge of the CDC. 
For electrons in the EC-EM we cut on the index of the 
most energetic tower, i^. Tower 15 covers 1.4 < 77 < 1.5 
with respect to the detector center and tower 25 covers 
2.4 < < 2.5. 



2. Quality Variables 

We test how well the shape of a cluster agrees with 
that expected for an electromagnetic shower by comput- 
ing a quality variable (x^) for all cell energies using a 
41-dimensional covariance matrix. The covariance ma- 
trix was determined from GEANT based simulations 

ii- 

To determine how well a track matches a cluster we ex- 
trapolate the track to the third EM layer in the calorime- 
ter and compute the distance between the extrapolated 
track and the cluster centroid in the azimuthal direction, 
As, and in the z-direction, Az. The variable 



As 

j7 



Az 



(20) 



quantifies the quality of the match. In the EC-EM z is 
replaced by r, the radial distance from the center of the 
detector. The parameters 5s — 0.25 cm, 5z = 2.1 cm, 
and Sr = 1.0 cm are the resolutions with which As, Az, 
and Ar are measured, as determined with the electrons 
from W ev decays. 

In the EC, electrons must have a matched track in the 
forward drift chamber. In the CC, we define "tight" and 
"loose" criteria. The tight criteria require a matched 
track in the CDC. The loose criteria do not require a 
matched track and help increase the electron finding ef- 
ficiency for Z ^ ee decays. 

The isolation fraction is defined as 



Er.. 



(21) 



where i?conc is the energy in a cone of radius R = 
•y/ A(/)2 -I- Arf = 0.4 around the direction of the elec- 
tron, summed over the entire depth of the calorimeter 
and Ecorc is the energy in a cone of R = 0.2, summed 
over the EM calorimeter only. 

Figure ^ shows the distributions of the three quality 
variables for electrons in the CC with the arrow showing 
the cut values. Table || summarizes the electron selection 
criteria. 



C. Electron Identification 



1. Fiducial Cuts 



D. Data Samples 



To ensure a uniform response we accept only electron 
candidates that are well separated in azimuth (A0) from 



The data were taken during the 1994-1995 Tevatron 
run. After the removal of runs in which parts of the 
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detector were not operating adequately, they amount to 
an integrated luminosity of about 82 pb~^. We select W 
decay candidates by requiring: 
Level 1: pp interaction 
Main Ring Veto 

EM trigger tower above 10 GeV 
Level 1.5: > 1 EM cluster above 15 GeV 
Level 2: electron candidate with pT > 20 GeV 

momentum imbalance > 15 GeV 
offline: > 1 tight electron candidate in CC 
prie) > 25 GeV 
Pt{i^) > 25 GeV 
UT < 15 GeV 
We select Z decay candidates by requiring: 
Level 1: pp interaction 

> 2 EM trigger towers above 7 GeV 
Level 1.5: > 1 EM cluster above 10 GeV 
Level 2: > 2 electron candidates with pT > 20 GeV 
offline: > 2 electron candidates 
PT{e) > 25 GeV 
70 < m(ee) < 110 GeV 
We accept Z ^ ee decays with at least one electron can- 
didate in the CC and the other in the CC or the EC. 
One CC candidate must pass the tight electron selec- 
tion criteria. If the other candidate is also in the CC 
it may pass only the loose criteria. We use the 2,179 
events with both electrons in the CC (CC/CC Z sam- 
ple) to calibrate the calorimeter response to electrons 
(Sec. VI). These events need not pass the Main Ring 
Veto cut because Main Ring background does not affect 
the EM calorimeter. The 2,341 events for which both 
electrons have tracks and which pass the Main Ring Veto 
(CC/CC-I-EC Z sample) serv e to calibrate the recoil mo- 
mentum response (Sec. VII). Table summarizes the 
data samples. 

Figure |lj shows the luminosity of the colliding beams 
during the W and Z data collection. 

On several occasions we use a sample of 295,000 ran- 
dom pp interaction events for calibration purposes. We 
collected these data concurrently with the W and Z sig- 
nal data, requiring only a pp interaction at Level 1. We 
refer to these data as "minimum bias events" . 



V. FAST MONTE CARLO MODEL 
A. Overview 

The fast Monte Carlo model consists of three parts. 
First we simulate the production of the W or Z boson by 
generating the boson four-momentum and other charac- 
teristics of the event like the ^-position of the interaction 
vertex and the luminosity. The event luminosity is re- 
quired for luminosity dependent parametrizations in the 
detector simulation. Then we simulate the decay of the 
boson. At this point we know the true px of the boson 
and the momenta of its decay products. We then apply 



a parameterized detector model to these momenta in or- 
der to simulate the observed transverse recoil momentum 
and the observed electron momenta. 



B. Vector Boson Production 

In order to specify completely the production dynam- 
ics of vector bosons in pp collisions we need to know the 
differential production cross section in mass Q, rapid- 
ity y, and transverse momentum qt of the produced W 
bosons. To speed up the event generation, we factorize 
this into 



dqj^dydQ dq^dy 



to generate Qt, Vi and Q of the bosons. 

For pp collisions, the vector boson production cross 
section is given by the parton cross section aij convo- 
luted with the parton distribution functions f[x,Q^) and 
summed over parton flavors 



d a 
dq^dy 



^ j dxi f dx2fi{xi,Q'^)fj{x2,Q'^) 

S3 



5{sxiX2 - Q ) 



,2- 
2\ Q 



dq^dy 



(23) 



Several authors i3C 



have computed ^j^^ 



us- 



ing a perturbative calculation 1 32 1 for the high-g^ regime 
and the CoUins-Soper resummation formalism [^,0 for 
the \ow-qT regime. We use the code provided by the 
authors of Ref. and the MRSA' parton distribution 
function s [p5| to compute the cross section. We evalu- 
ate Eq. |23| separately for interactions involving at least 
one valence quark and for interactions involving two sea 
quarks. 

The parton cross section is given by 



d^a 
dq^dy 



A'KS 



d%e'^^-^ ■ W{h) xe-^ + Y 



(24) 



where ctq is the tree-level cross section, s is the parton 
center-of-mass energy, and b is the impact parameter in 
transverse momentum space. W and Y are perturbative 
terms and S parameterizes the non-perturbative physics. 
In the notation of Ref. 1 30 



S = 



91 



'''""[Wo 



gigs ln(100a;iX2)6 



(25) 



where Qo is a cut-off parameter, xi and X2 are the mo- 
mentum fractions of the initial state partons. The param- 
eters Qi, ff2 , and gs have to be determined experimentally 

(Seclvn^). 
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We use a Breit-Wigner curve with mass dependent 
width for the hne shape of the W boson. The intrin- 
sic width of the W isrw = 2.062 ±0.059 GeV ||]. The 
hne shape is skewed due to the momentum distribution 
of the quarks inside the proton and antiproton. The mass 
spectrum is given by 



(26) 



Ml 



We call 

the parton luminosity. To evaluate it we generate W — > 
ev events using the herwig Monte Carlo event generator 
p7t , interfaced with pdflib Q, and select the events 
subject to the same kinematic and fiducial cuts as for 
the W and Z samples with all electrons in CC. We plot 
the mass spectrum divided by the intrinsic line shape of 
the W boson. The result is proportional to the parton 
luminosity and we parameterize the spectrum with the 
function [ITl 



-PQ 



(28) 



Table [II shows f3 for W and Z events for some modern 
parton distribution functions. The value of P depends 
on the rapidity distribution of the W bosons, which is re- 
stricted by the kinematic and fiducial cuts that we impose 
on the decay leptons. The values of /? given in Table 111 
are for the rapidity distributions of W and Z bosons that 



satisfy the kinematic and fiducial cuts given in Sec. IV 



The uncertainty in /3 is about 0.001, due to Monte Carlo 
statistics and uncertainties in the acceptance. 

To generate the boson four-momenta we treat da/dQ 
and d^a/dq^dy as probability density functions and pick 
Q from the former and a pair of y and qt values from 
the latter. For a fraction fss we use d^a/dq^dy for inter- 
actions between two sea quarks. Their helicity is +1 or 
— 1 with equal probability. For the remaining W bosons 
we use d^a/dq^dy for interactions involving at least one 
valence quark. They always have helicity —1. Finally, 
we pick the ^-position of the interaction vertex from a 
Gaussian distribution centered at z = with a standard 
deviation of 25 cm and a luminosity for each event from 
the histogram in Fig. ^ 



C. Vector Boson Decay 

At lowest order the W boson is fully polarized along 
the beam direction due to the V-A coupling of the 
charged current. The resulting angular distribution of 
the charged lepton in the W rest frame is given by 



da 



where A is the helicity of the W with respect to the proton 
direction, q is the charge of the lepton, and 9* is the angle 
between the charged lepton and proton beam directions 
in the W rest frame. The spin of the W points along the 
direction of the incoming antiquark. Most of the time 
the quark comes from the proton and the antiquark from 
the antiproton, so that A = — 1. Only if both quark and 
antiquark come from the sea of the proton and antiproton 
is there a 50% chance that the quark comes from the 
antiproton and the antiquark from the proton and in that 
case A = 1 (Fig. |l^). We determine the fraction of sea- 
sea interactions, /ss, using the parameterizations of the 
parton distribution functions given in pdflib |Q . 

When 0{as) processes are included, the boson acquires 
finite pT and Eq. E9| is changed to 



do- 



d cos ^cs 



cx (l + ai{qT) coaOcs + ct2{qT) coa^ 9cs) 



(30) 



d cos 9* 



(X (1 - Agcosr)^ 



(29) 



for bosons with A = — 1 and after integration over 
(j>. The angle ^cg in Eq. 30 is now defined in the CoUins- 
Soper frame p0[ | . The values of ai and a2 as a function 
of transverse boson momentum have been calculated at 
0{a'^) 1^ and are shown in Fig. |l^. We have imple- 
mented the angular distribution given in Eq. ^ in the 
fast Monte Carlo. The effect is smaller if the W bosons 
are selected with ut < 15 GeV than for ut < 30 GeV. 
The angular distribution of the leptons from Z —^ ee de- 
cays is also generated according to E g. po[ but with ai 
and a2 computed for Z ^ ee decays [p9|. 

To check whether neglecting the correlations between 
the mass and the other parameters in Eq. ^ introduces 
an uncertainty, we use the herwig program to generate 
W ^ ev decays including the correlations neglected in 
our model. We apply our parameterized detector model 
to them and fit them with probability density functions 
that were generated without the correlations. The fitted 
W mass values agree with the W mass used in the Monte 
Carlo generation within the statistical uncertainties of 25 
MeV. 

Radiation from the decay electron or the W boson bi- 
ases the mass measurement. If the decay electron radi- 
ates a photon and the photon is well enough separated 
from the electron so that its energy is not included in 
the electron energy, or if an on-shell W boson radiates 
a photon and therefore is off-shell when it decays, the 
measured mass is biased low. We use the calculation of 
Ref. |^l|] to generate W ev^ decays. The calculation 
gives the fraction of events in which a photon with energy 
i?(7) > Eq is radiated, and the angular distribution and 
energy spectrum of the photons. Only radiation from the 
decay electron and the W boson, if the final state W is 
off-shell, is included to order a. Radiation by the initial 
quarks or the W , if the final W is on-shell, does not affect 
the mass of the ev pair from the W decay. We use a min- 
imum photon energy Eq = 50 MeV, which means that in 
30.6% of all W decays a photon with Eij) > 50 MeV 
is radiated. Most of these photons are emitted close to 
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the electron direction and cannot be separated from the 
electron in the calorimeter. For Z — > ee decays there is 
a 66% probabihty that any one of the electrons radiates 
a photon with Ei-y) > 50 MeV. 

The separation of the electron and photon in the lab 
frame is 



ARiej) = V(0(e) - </>(7))2 + ('7(e) - '?(7))= 



(31) 



Figure ^ shows the calculated distribution of photons as 
a function of Ai?(e7). The shaded histogram in the figure 
shows the photons that are reconstructed as separate ob- 
jects. If the photon and electron are close together they 
cannot be separated in the calorimeter. The momen- 
tum of a photon with Ai?(e7) < Rq is therefore added 
to the electron momentum, while for Ai?(e7) > Rq & 
photon is considered separated from the electron and its 
momentum is added to the recoil momentum. We use 
Rq = 0.3, which is the approximate size of the window in 
which the electron energy is measured. This procedure 
has been verified to give the same results as an explicit 
GEANT simulation of radiative W decays. In only about 
3.5% of the W ^ ev decays does the photon separate far 
enough from the electron, i.e. Ai?(e7) > i?Oi to cause a 
mismeasurement of the transverse mass. 

W boson decays through the channel W ^ tv ^ evVv 
are topologically indistinguishable from W ev de- 
cays. We therefore include these decays in the W de- 
cay model, properly accounting for the polarization of 
the tau leptons in the decay angular distributions. The 
fraction of W bosons that decay in this way is B{t 
evV)/ (1 -f B{t euV)) = 0.151. 

We let the generated W bosons decay with an angular 
distribution corresponding to their helicity. For 15.1% of 
the W bosons the decay is to tv evvv. For 30.6% 
of the remaining W bosons a photon is radiated. For 
66% of the Z bosons the decay is to e+e"7 and for the 
remainder to e^e~ . 



D. Detector Model 



The transverse energy measurement depends on the 
measurement of the electron direction as well. We de- 
termine the shower centroid position by intersecting the 
line defined by the event vertex and the electron direction 
with a cylinder coaxial with the beam and 91.6 cm in ra- 
dius (the radial center of the EM3 layer). We then smear 
the azimuthal and z-coordinates of the intersection point 
by their resolutions. We determine the z-coordinate of 
the center of gravity of the CDC track by intersecting the 
same line with a cylinder of 62 cm radius, the mean radial 
position of all delay lines in the CDC, and smearing by 
the resolution. The measured angles are then obtained 



from the smeared points as described in Section IV B 1 



The model for the particles recoiling against the W 
has two components: a "hard" component that models 
the pt of the W , and a "soft" component that mod- 
els detector noise and pile-up. Pile-up refers to the ef- 
fects of additional pp interactions in the same or previous 
beam crossings. For the soft component we use the trans- 
verse momentum balance jfx from a minimum bias event 
recorded in the detector. The observed recoil pt is then 
given by 



UT = -(RrccgT + 0-rcc ' X)qT 

-Au||(/:,M||)pT(e) 

+ amb/T, 



(33) 



where qT is the generated value of the boson transverse 
momentum, Rj-oc is the (in general momentum depen- 
dent) response, tXrec is the resolution of the calorimeter, 
Au|| is the transverse energy flow into the electron win- 
dow (parameterized as a function of luminosity C and 
U||), and ttnib is a correction factor that allows us to ad- 
just the resolution to the data. The quantity Aii|| is 
different from the energy added to the electron, Ai?, be- 
cause of the zero-suppression in the calorimeter readout. 

We simulate selection biases due to the trigger require- 
ments and the electron isolation by accepting events with 
the estimated efficiencies. Finally, we compute all the de- 
rived quantities from these observables and apply fiducial 
and kinematic cuts. 



The detector simulation uses a parameterized model 
for response and resolution to obtain a prediction for the 
distribution of the observed electron and recoil momenta. 

When simulating the detector response to an electron 
of energy Eq, we compute the observed electron energy 
as 

E{e)^a^MEfi + AE{C,u\\) + G^M-X, (32) 

where q;em is the response of the electromagnetic 
calorimeter, /S.E is the energy due to particles from the 
underlying event within the electron window (parameter- 
ized as a function of luminosity C and U||), ctem is the 
energy resolution of the electromagnetic calorimeter, and 
X is a random variable from a normal parent distribution 
with zero mean and unit width. 



VI. ELECTRON MEASUREMENT 

A. Angular Resolutions 

The resolution for the z-coordinate of the track center 
of gravity, ztrk, is determined from the Z — > ee sample. 
Both electrons originate from the same interaction ver- 
tex and therefore the difference between the interaction 
vertices reconstructed from the two electrons separately, 
Zvt^{ei) — Zvtx(e2), is a measure of the resolution with 
which the electrons point back to the vertex. The points 
in Fig. |2^ show the distribution of Zvtx(ei) — Zvtx(e2) ob- 
served in the CC/CC Z sample with tracks required for 
both electrons. 
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A Monte Carlo study based on single electrons gener- 
ated with a GEANT simulation shows that the resolution 
of the shower centroid algorithm can be parameterized 

as 

fT(2cai) = (a + 5A(e)) + (c + dA(e))zeai, (34) 

where A(e) = \0{e) - 90°|, a = 0.33 cm, b = 5.2 x 10~^ 
cm, c = 4.2 X 10""', and d = 7.5 x 10"^. We then tune 
the resolution function for ztrk in the fast Monte Carlo 
so that it reproduces the shape of the Zvtx(ei) — z^tx{^2) 
distribution observed in the data. We find that a reso- 
lution function consisting of two Gaussians 0.31 cm and 
1.56 cm wide, with 6% of the area under the wider Gaus- 
sian, fits the data well. The histogram in Fig. |2l| shows 
the Monte Carlo prediction for the best fit, normalized 
to the same number of events as the data. The W mass 
measurement is very insensitive to these resolutions. The 
uncertainties in the resolution parameters cause less than 
5 MeV uncertainty in the fitted W mass. 

The calibration of the z-position measurements from 
the CDC and calorimeter is described in Appendix 
We quantify the calibration uncertainty in terms of scale 
factors acDC = 0.988 ±0.001 and acc = 0.9980 ± 0.0005 
for the z-coordinate. The uncertainties in these scale 
factors lead to a finite uncertainty in the W mass mea- 
surement. 



B. Underlying Event Energy 

The energy in an array of 5x5 towers in the four EM 
layers and the first FH layer around the most energetic 
tower of an electron cluster is assigned to the electron. 
This array contains the entire energy deposited by the 
electron shower plus some energy from other particles. 
The energy in the window is excluded from the compu- 
tation of ut. This causes a bias in the component of 
ut along the direction of the electron. For pt{W) ^ AIw 



niT ~ 2pT(e) + u 



(35) 



so that this bias propagates directly into a bias in the 
transverse mass. We call this bias Aw||. It is equal to 
the momentum flow observed in the EM and first FH 
sections of a 5x5 array of calorimeter towers. 

We use the W and Z data samples to measure Au||. 
For every electron in the W and Z samples we compute 
the energy flow into an azimuthal ring of calorimeter tow- 
ers, 5 towers wide in 77 and centered on the tower with 
the largest fraction of the electron energy. For every elec- 
tron we plot the transverse energy flow into one-tower- 
wide azimuthal segments of this ring as a function of the 
azimuthal separation |A0| between the center of the seg- 
ment and the electron shower centroid. The energy flow 
^ £^1x5 is computed as the sum of all energy deposits in 
the four EM layers and the first FH layer in the 1x5 tower 
segment. Figure shows the transverse energy flow 



X] -''^1x5/ cosh 77(e) versus |A(/)| for the electrons in the 
W sample with ut < 15 GeV. For small \A(j)\ we see the 
substantial energy flow from the electron shower and for 
larger \A(j)\ the constant noise level. The electron shower 
is contained in a window of |A0| < 0.2. We estimate the 
energy flow into the 5x5 tower window around the elec- 
tron from the energy flow into segments of the azimuthal 
ring with \A(j)\ > 0.2. The level of energy flow is sensitive 
to the isolation cut. The region 0.2 < |A<^| < 0.4, which 
is used for the isolation variable, is maximally biased by 
the cut; the region, 0.4 < |A0| < 0.6, which is close to 
the electron but outside the isolation region, is minimally 
biased. We expect the energy flow under the electron to 
lie somewhere in between the energy flow into these two 
regions. We therefore compute Au|| based on the aver- 
age transverse energy flow into both regions and assign a 
systematic error equal to half the difference between the 
two regions. We repeat the same analysis for the elec- 
trons in the CC/CC Z sample. The results are tabulated 
in Table 1^. We find Awy = 479 ±2{stat)±6{syst) MeV 
for W events with ut < 15 GeV. For the Z sample Am|| 
is 11 ± 7 MeV lower. Figure p3| shows the spectrum of 



Aii 



At higher luminosity the average number of interac- 
tions per event increases and therefore Au|| increases. 
This is shown in Fig. The mean value of Auy in- 
creases by 11.2 MeV per 10^"cm~2s~^ The underlying 
event energy flow into the electron window also depends 
on uy. Figure p5| shows (Au||(0,M||)), the mean value for 
Am II corrected back to zero luminosity, as a function of 
In the fast Monte Carlo model a value Aity is picked 
from the distribution shown in Fig. ^ for every event 
and then corrected for wy and luminosity dependences. 

The measured electron energy is biased upwards by the 
additional energy AE in the window from the underly- 
ing event. AE is not equal to Au|| because the additional 
energy deposited by the electron may lift some cells that 
would have been zero-suppressed in the calorimeter read- 
out above the zero-suppression threshold. Therefore 



AE = A-itii 



(36) 



where Apod — 212 ± 25 MeV is a correction for the 
pedestal shift introduced by the zero-suppression in the 
calorimeter readout. This is determined by superimpos- 
ing single electrons simulated with a GEANT simulation 
on minimum bias events that were recorded without zero- 
suppression in the calorimeter readout. Most of this bias 
cancels in the W to Z mass ratio so that the W mass 
measurement is not sensitive to Aped- 



C. M|[ Efficiency 

The efficiency for electron identification depends on 
their environment. Well- isolated electrons are identified 
correctly more often than electrons near other particles. 
Therefore W decays in which the electron is emitted in 
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the same direction as the particles recoiling against the 
W are selected less often than W decays in which the 
electron is emitted in the direction opposite the recoiling 
particles. This causes a bias in the lepton px distribu- 
tions, shifting pt{s.) to larger values and pt{i^) to lower 
values, whereas the rriT distribution is only slightly af- 
fected. 

We estimate the electron finding efficiency as a func- 
tion of it|| by superimposing Monte Carlo electrons, sim- 
ulated using the GEANT program, onto the events from 
our W signal sample. We use the W sample in order 
to ensure that the underlying event is correctly modeled. 
The sample of superimposed electrons, which are spa- 
tially separated from the electron that is already in the 
event, matches the data well. It is important that the 
superimposed sample model the transverse shower shape 
and isolation well, because these are the dominant ef- 
fects that cause the efficiency to vary with u\\. Figure |26| 
shows the transverse shower profile of the superimposed 
electron sample and the electron sample from W decays. 
Figure ^ shows the distribution of the isolation for the 
two electron samples in five U|| regions. Figure com- 
pares the mean isolation versus U|| for the two samples. 

We then apply the shower shape and isolation cuts used 
to select the W signal sample and determine the fraction 
of the electrons in the superimposed samples that pass 
all requirements as a function of This efficiency is 
shown in Fig. |2^. The line is a fit to a function of the 
form 



tron energy 



e(w|l 



1 for U|| < itQ 

1 — s(uj| — uo) otherwise. 



(37) 



The parameter Eq is an overall efficiency which is inconse- 
quential for the W mass measurement, uq is the value of 
M|| at which the efficiency starts to decrease as a function 
of M|| , and s is the rate of decrease. We obtain the best fit 
for Mo = 3.85 ± 0.55 GeV and s = 0.013 ± 0.001 GeV-i. 
These two values are strongly correlated. The errors ac- 
count for the finite number of superimposed Monte Carlo 
electrons. 



D. Electron Energy Response 

Equation |7 relates the reconstructed electron energy 
to the recorded calorimeter signals. Since the values 
for the constants were determined in a different setup, 
we determine the offset (5em and a scale asM, which es- 
sentially modifies A, in situ with collider data for reso- 
nances that decay to electromagnetically showering par- 
ticles: 7r° 77, J jil) ^ e+e", and Z ee. We use 7r° 
and J I ip signals from an integrated luminosity of approx- 
imately 150 nb~^, accumulated during dedicated runs 
with low pt thresholds for EM clusters in the trigger. 

The fast Monte Carlo predicts the reconstructed elec- 



-'EM 



(38) 



where Eq is the generated electron energy. To determine 
(5em and aEM, we compare the observed resonances and 
Monte Carlo predictions as a function of a^M and Sem- 
The photons from the decay of 7r°s with pT > 1 GeV 
cannot be separated in the calorimeter. There is about 
a 10% probability for each photon to convert to an e+e~ 
pair in the material in front of the CDC. If both photons 
convert we can identify tt" decays as EM clusters in the 
calorimeter with two doubly-ionizing tracks in the CDC. 
We measure the 7r° energy E{Tr°) in the calorimeter and 
the opening angle lu between the two photons using the 
two tracks. This allows us to compute the "symmetric 
mass" 



''sym 



1 



cos w 



(39) 



which is equal to the invariant mass if both photons have 
the same energy, and is larger for asymmetric decays. 
Figure 30 shows the background subtracted spectrum of 



''sym 



for 7r° candidates in the CC-EM superimposed with 
a Monte Carlo prediction of the line shape. 

Figure ^ shows the invariant mass spectrum of dielec- 
tron pairs in the J/ip mass region. The smooth curve is 
the fit to a Gaussian line shape above the background 
predicted using a sample of EM clusters without CDC 
tracks. After correction for underlying event effects we 
measure a mass of 3.03 ± OM{stat)_± 0.19{syst) GeV. A 
Monte Carlo simulation of pp ^ bb + X, b J/ip + X 
tells us that we expect to observe a mass 



mobs = aEMmj/^ + 0.56 Sem- 



(40) 



Together with our measurement of Wobs) this restricts 
the allowed parameter space for oem and Sem- The 7r° 
and J/ip analyses are described in detail in Ref. 
Figure |3^ shows the 68% confidence level contours in q;em 
and Sem obtained from these data. 

Fixing the observed Z boson mass to the measured 
value (Eq. ||) correlates the values allowed for uem and 
Sem- For a given Sem we determine uem so that the po- 
sition of the Z peak predicted by the fast Monte Carlo 
agrees with the data. To determine the scale factor 
that best fits the data, we perform a maximum likeli- 
hood fit to the m(ee) spectrum between 70 GeV and 
110 GeV. In the resolution function we allow for an 
exponential background shape whose slope is fixed to 
-0.037 ± 0.002 GeV~\ the value obtained from a sam- 
ple of events with two EM clusters that fail the electron 
quality cuts (Fig. |3^). The background normalization is 
allowed to float in the fit. This is sufficient, together with 
the tt'^ and J/ip data, to determine both uem and Sem- 
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Without relying on the low energy data at all we can 
extract aEM and Jem from the Z data alone. The elec- 
trons from Z decays are not monochromatic and there- 
fore we can make use of their energy spread to constrain 
Q!EM and (5em simultaneously. For (5em ^ -^'(ei) + ^{^"i) 
we can write 



m(ee) = qiemMz + /zfeM 



(41) 



where ]z — (i?(ei)-|-i?(e2))(l — cosLi;)/m(ee) and uj is the 
opening angle between the two electrons. We plot m{ee) 
versus fz (Fig. 33 ) and compare it with the Monte Carlo 
predictions for the allowed values of a^M and Sem using 
a binned maximum likelihood fit. 

Using the constraints on oem and Sem from the Z data 
alone we obtain the contour labeled "Z" in Fig. |3j and 
Sem = 0.02 ±0.36 GeV. The uncertainty in this measure- 
ment of Sem is dominated by the statistical uncertainty 
due to the finite size of the Z sample. 

The combined constraint from all three resonances is 
shown by the thick contour in Fig. |4[ The 7r° and J/tp 
contours essentially fix Sem, independent of ckem- The re- 
quirement that the Z peak position agree with the known 
Z boson mass correlates q;em and iJem. The contours in 
Fig. ^ reflect only statistical uncertainties. The uncer- 
tainty in the tt^ and J/ip contours is dominated by sys- 
tematic effects in the underlying event corrections and 
the deviation of the test beam data from the assumed 
response at low energies. The double arrow in Fig. ^ 
represents the systematic uncertainty in Sem- We deter- 
mine 



Sem = -0.161°:^? 



(42) 



Figure ^ shows the m(ee) spectrum for the CC/CC 
Z sample and the Monte Carlo spectrum that best fits 
the data for Sem = —0.16 GeV. The for the best fit 
to the CC/CC m[ee) spectrum is 33.5 for 39 degrees of 
freedom. For 



aEM = 0.9533 ±0.0008 



(43) 



the Z peak position is consistent with the known Z boson 
mass. The error reflects the statistical uncertainty and 
the uncertainty in the background normalization. The 
background slope has no measurable effect on the result. 

If we split the CC/CC Z sample into events with two 
tight electrons and events with a tight and a loose elec- 
tron and fit them separately using the value of uem given 
in Eq. ^ we obtain 

Mz = 91.206 ± 0.086 GeV (tight/tight sample); (44) 
Mz = 91.145 ± 0.148 GeV (tight/loose sample). (45) 

Figures |3^ (a) and (b) show the corresponding spectra 
and fits. 



E. Electron Energy Resolution 

Equation ^ gives the functional form of the electron 
energy resolution. We take the intrinsic resolution of 
the calorimeter, which is given by the sampling term 
Sem, from the test beam measurements. The noise term 
riEM is represented by the width of the Ai? distribution 
(Fig. ^3|). We measure the constant term cem from the 
Z line shape of the data. We fit a Breit-Wigner con- 
voluted with a Gaussian, whose width characterizes the 
dielectron mass resolution, to the Z peak. Figure ^ 
shows the width (Tm(ee) of the Gaussian fitted to the Z 
peak predicted by the fast Monte Carlo as a function of 
Cem- The horizontal lines indicate the width of the Gaus- 
sian fitted to the CC/CC Z sample and its uncertainties, 
1.75 ± 0.08 GeV. We find that Monte Carlo and data 
agree if cem = 0.0115lg'QQ3g, as indicated by the arrows 
in Fig. The measured Z mass does not depend on 
Cem- 



VII. RECOIL MEASUREMENT 
A. Recoil Momentum Response 

The detector response and resolution for particles re- 
coiling against a W boson should be the same as for par- 
ticles recoiling against a Z boson. For Z — s- ee events, 
we can measure the transverse momentum of the Z from 
the e"'"e~ pair, pT(ee), into which it decays and from the 
recoil momentum ut in the same way as for W ev 
events. By comparing pT{ee) and ut we calibrate the 
recoil response relative to the electron response. 

The recoil momentum is carried by many particles, 
mostly hadrons, with a wide momentum spectrum. Since 
the response of calorimeters to hadrons tends to be non- 
linear and the recoil particles are distributed all over the 
calorimeter, including module boundaries with reduced 
response, we expect a momentum dependent response 
function with values below unity. In order to fix the func- 
tional form of the recoil momentum response, we study 
the response predicted by a Monte Carlo Z ee sample 
obtained using the herwig program and a GEANT-based 
detector simulation. We project the reconstructed trans- 
verse recoil momentum onto the direction of motion of 
the Z and define the response as 



f^rcc 



\UT ■ qr] 
|9t| 



(46) 



where qx is the generated transverse momentum of the 
Z boson. Figure |3^ shows this response as a function of 
qt- a response function of the form 



R-rcc = ttroc + A'cc log (^t/GbV) 



(47) 



fits the response predicted by GEANT with a^^c = 0.713 ± 
0.006 and /3rec = 0.046 ±0.002. This functional form also 
describes the jet energy response of the D0 calorimeter. 
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To measure the recoil response from the collider data 
we use the CC/CC+EC Z sample. We allow one of the 
leptons from the Z — > ee decay to be in the CC or the 
EC, so that the rapidity distribution of the Z bosons 
approximates that of the W bosons. We require both 
leptons to satisfy the tight electron criteria. This reduces 
the background for the topology with one lepton in the 
EC. We also require the Main Ring Veto as for the W 
sample (Sec. iV). 

We project the transverse momenta of the recoil, u^, 
and the Z as measured by the two electrons, pxiee), on 
the inner bisector of the electron directions (77-axis), as 
shown in Fig. [T4| By projecting the momenta on an axis 
that is independent of any energy measurement, noise 
contributions to the momenta average to zero and do not 
bias the result. We bin the data in p^(ee) and plot the 



mean of the sum of the two projections. 



versus the mean of p,,(ee) (Fig. p^). We perform a two- 
dimensional fit for the two parameters by compar- 
ing the data to predictions of the fast Monte Carlo for 
different values of a^cc a-nd /3ioc- Figure 39 also shows 
the prediction of the Monte Carlo for the values of the 
parameters that give the best fit. Figure ^ shows the 
contour for = Xo + 1- The best fit (xo = 5 for 8 de- 
grees of freedom) is achieved for aioc — 0.693 ±0.060 and 
/3roc = 0.040 ± 0.021. The two parameters are strongly 
correlated with a correlation coefficient p = —0.979. 



B. Recoil Momentum Resolution 

We parameterize the resolution for the hard compo- 
nent of the recoil as 



UT, 



(48) 



where Srec is a tunable parameter. 

The soft component of the recoil is modeled by the 
transverse momentum balance j^t from minimum bias 
events, multiplied by a correction factor amb (Eq. 33). 
This automatically models the effects of detector reso- 
lution and pile-up. To model the pile-up correctly as 
a function of luminosity, we need to take the minimum 
bias events at the same luminosity as the W events. At 
a given luminosity the mean number of interactions in 
minimum bias events is always smaller than the mean 
number of interactions in W events. To model the detec- 
tor resolution correctly, the minimum bias events must 
have the same interaction multiplicity spectrum as the 
W events. We therefore weight the minimum bias events 
so that their interaction multiplicity approximates that 
of the W events. As a measure of the interaction multi- 
plicity on an event-by-event basis, we use the multiplicity 
of vertices reconstructed from the tracks in the CDC and 
the timing structure of the Level hodoscope signals Q . 

We tune the two parameters Srcc and amb using the 
CC/CC-HEC Z sample. The width of the spectrum of 



the Ty-balance, w^/Rrec +Pri{s&), is a measure of the re- 
coil momentum resolution. Figure ^ shows this width cr^ 
as a function of p,,(ee). The contribution of the electron 
momentum resolution to the width of the 77-balance is 
negligibly small. The contribution of the recoil momen- 
tum resolution grows with p^(ee) while the contribution 
from the minimum bias j^t is independent of (ee). This 
allows us to determine s^cc and Ofmb simultaneously and 
without sensitivity to the electron resolution by compar- 
ing the width of the jy-balance predicted by the Monte 
Carlo model with that observed in the data in bins of 
Prf{ee). We perform a fit comparing Monte Carlo and 
collider data. Figure ^ shows contours of constant in 
the ainb-Srcc plane. The best agreement (xq = 10.3 for 8 
degrees of freedom) occurs for Srec = 0.49 ± 0.14 GeV^^^ 
and ttnib = 1.032 ± 0.028 with a correlation coefficient 
p — —0.60 for the two parameters. The ^-balance, 
u^/Ricc +P5(ee), is more sensitive to the electron mo- 
mentum resolution and is affected by changes in s^cc and 
amb in the same way. We use it as a cross check only. 

Figure ^ shows the spectrum of u,,/Ri-cc +Pj;(ee) from 
the CC/CC-I-EC Z data sample and from the fast Monte 
Carlo with the tuned recoil resolution and response pa- 
rameters. Figure ^ shows the corresponding distribu- 
tions for u^/Ricc +Pc(ee). In both cases the agreement 
between data and Monte Carlo simulation is good. A 
Kolmogorov-Smirnov test |Q gives confidence levels of 
K — 0.33 and 0.37 that the Monte Carlo and data spectra 
derive from the same parent distribution. A test gives 
= 25 and 37, respectively, for 40 bins. 

Figure ^ shows the overall energy flow transverse 
to the beam direction measured by the sum Sx — 
J2i Ei sin Oi over all calorimeter cells except cells belong- 
ing to an electron cluster. For W events {St) = 98.7±0.3 
GeV and for Z events {St) = 91.0 ± 0.9 GeV. Increased 
transverse energy flow leads to a worse recoil momentum 
resolution and therefore we need to correct the value of 
amb for the W sample to account for this difference. Fig- 
ure ^ relates transverse energy flow St to resolution ctt 
for a minimum bias event sample. The resolution for 
measuring transverse momentum balance along any di- 
rection is 

ctt(S't) = 1.42 GeV + 0.15\/5^^GeV -f 0.007S't (49) 

for minimum bias events. The different energy flows 
in W and Z events lead to a correction to amb of 
crT(98.7 GeV)/crT(91.0 GeV) = 1.03 ± 0.01. The un- 
certainty reflects the uncertainties in the determination 
of {St)- This uncertainty does not correlate with Srcc- 

Z bosons are not intrinsically produced with less en- 
ergy flow in the underlying event than W bosons. Rather, 
the requirement of two reconstructed isolated electrons 
biases the event selection in the Z sample towards events 
with lower energy flow compared to the events in the W 
sample which have only one electron. We demonstrate 
this by loosening the electron identiflcation requirements 
for one of the electrons in the Z sample. We use events 
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that were collected using less restrictive trigger condi- 
tions for which at Level 2 only one of the electron candi- 
dates must satisfy the shape and isolation requirements. 
We find that if all electron quality cuts are removed for 
one electron St increases by 7%, consistent with the ratio 
of the St values in the W and Z samples. 

C. Comparison with W Data 



g, = 0.1lt°:°^ GeV^; 

32 = 0.58i|^:2 GeV^ (50) 

for mass cut-off Qo = 1-6 GcV in Eq. || and CTEQ2M 
parton distribution functions, by fitting Drell-Yan and Z 
data at different values of Q^. We further constrain these 
parameters using our much larger Z data sample. 



B. Determination of g2 from Z — > ee Data 



'I 

—0.64 ±0.03 GeV and for the Monte Carlo prediction in- 
cluding backgrounds it is —0.61 ± 0.01 GeV, in very good 
agreement. This is important because a bias in W|| would 
translate into a bias in the determination of rriT (Eq. 35 ) . 
The agreement means that recoil momentum response 
and resolution and the M|| efficiency parameterization de- 
scribe the data well. Figures ^-[50| show u±, ut, and the 
azimuthal difference between electron and recoil direc- 
tions from Monte Carlo and W data. The Kolmogorov- 
Smirnov probabilities for Figs. |47-50 
0.16, and 0.11, respectively. 



are k = 0.15, 0.38, 



We compare the recoil momentum distribution in the 
W data to the predictions of the fast Monte Carlo, which 
includes the parameters determined in this section and 

Sec. 0. Figure ^ compares the u,, spectra from Monte The pt of Z bosons can be measured more precisely 
Carfo and W data. The mean m for the W data is ^^an the pt of W bosrais by using the e+e pairs from 

their decays. Figure pll shows the pT{ee) spectrum ob- 
served in the data. 

To reduce the background contamination of the sam- 
ple, the invariant mass of Z candidates must be within 
10.5 GeV of the Z peak position. This mass window re- 
quirement reduces the background fraction to 2.5%, as 
determined from the dielectron invariant mass spectrum. 
As such it includes a contribution from Drell-Yan e~^e~ 
production, which has a pT spectrum similar to the sig- 
nal and should not be counted as background in this case. 
To account for this uncertainty we assign an error to the 
background fraction of ±2.5%. 

The shape of the background is fixed by a sample of 
events with two electromagnetic clusters which pass the 
same kinematic requirements as our Z ee sample, but 
fail the electron identification cuts (sample 1). As a 
cross-check we also use events with two jets, each with 
more than 70% of its energy in the EM calorimeter (sam- 
ple 2). Parameterizations of the two background shapes 
are shown in Fig. Their difference is taken to be the 
uncertainty in background shape. 

We use the fast Monte Carlo model to predict the 
PT{ee) spectrum from Z — )■ ee decays for different sets 
of parameter values. The fast Monte Carlo simulates 
the detector acceptance and resolution as discussed in 
the previous sections. Figure |5^ shows the pT{ee) spec- 
tra predicted by the fast Monte Carlo for MRSA' parton 
distribution functions and three values of g2, with gi and 
53 fixed at the values given in Eq. |5^. 

The dominant effect of varying 52 is to change the 
mean boson p-p. Properly normalized and with the back- 
ground contribution added, we use these distributions 
as probability density functions to perform a maximum 
likelihood fit for g2. For a set of discrete values of 52 
we compute the joint likehhood L of the observed pr(ee) 
spectrum. We then fit logi as a function of 52 with a 
third order polynomial. The maximum of the polynomial 
gives the fitted value of 52- The value of 52 has to be fit in- 
dependently for each parton distribution function choice. 
We perform fits for four choices of parton distribution 
functions: MRSA', MRSD-', CTEQ2M, and CTEQ3M. 
We fit the spectrum over the range pT(ee) < 15 GeV, 



VIII. CONSTRAINTS ON THE W BOSON pr 
SPECTRUM 

A. Parameters 

Since we cannot reconstruct a Lorentz invariant mass 
for W ev decays, knowledge of the transverse momen- 
tum distribution of the W bosons is necessary to measure 
the mass from the kinematic distributions. Theoretical 
calculations provide a formalism to describe the boson pT 
spectrum, but it includes phcnomcnological parameters 
gi, g2, and 03, which need to be determined experimen- 
tally (Sec. VB). In addition, the boson spectrum also 
depends on the choice of parton distribution functions 
and Aqcd- 

We can measure the W boson pT spectrum only indi- 
rectly by measuring ut, the pt of all particles that recoil 
against the W boson. Momentum conservation requires 
the W boson pT to be equal and opposite to ut- The pre- 
cision of the Ut measurement is insufficient, especially for 
small Ut, to constrain the W spectrum as tightly as is 
necessary for a precise W mass measurement. 

We therefore have to find other data sets to constrain 
the model. The formalism that describes the pT spec- 
trum of the W bosons has to simultaneously describe the 
Pt spectrum of Z bosons and the dilepton pT spectrum 
from Drell-Yan production with the same model param- 
eter values. The authors of Ref. ISOll find 
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which corresponds to the range accepted by the W selec- 
tion cuts. The fits describe the data well. Table lists 
the fitted values for 52 for the different parton distribu- 
tion function choices. The result of the CTEQ2M fit is 
in good agreement with the value in Eq. |50| . 

We estimate systematic uncertainties in the 32 fit by 
running the fast Monte Carlo with different parameter 
values and refitting the predicted priee) spectrum with 
the nominal probability density functions. The uncer- 
tainties in electron momentum response and resolution, 
M|| efficiency parametrization, fiducial cuts, model of ra- 
diative decays, and background translate into a system- 
atic uncertainty in g2 of 0.05 GeV^. 

As a cross-check we also fit the spectrum of the az- 
imuthal separation A0(ee) between the two electrons to 
constrain (72- The A(/)(ee) spectrum has smaller system- 
atic uncertainties but less statistical sensitivity to g2 than 
the A(j){ee) spectrum. In Table ^ we also quote the re- 
sults for (72 from a fit to the A0(ee) spectrum. 

The Monte Carlo prediction for the fitted 52 value us- 
ing MRSA' parton distribution functions is superimposed 
as a smooth curve on Fig. The Kolmogorov-Smirnov 
probability that the two distributions are from the same 
parent distribution is k = 0.72 and the is 25.5 for 29 
degrees of freedom. Both of these tests indicate a good 
fit. We use this model to compute the probability density 
functions for the final fits to the kinematic spectra from 
the W sample. 



IX. BACKGROUNDS 

A. W —> TV —> evVv 

The decay W ^ tv ^ evvv is topologically indistin- 
guishable from W ev. It is included in the fast Monte 
Carlo simulation (Sec. 0). This decay is suppresse d by 
the branching fraction for r evv, (17.83 ± 0.08)% @, 
and by the lepton -pT cuts. It accounts for 1.6% of events 
in the W sample. 



B. Hadronic Background 

QCD processes can fake the signature oi a W ^ ev 
decay if a hadronic jet fakes the electron signature and 
the transverse momentum balance is mismeasured. 

We estimate this background from the '^t spectrum of 
events with an electromagnetic cluster. Electromagnetic 
clusters in events with low j^t are almost all due to jets. 
A fraction satisfy our electron selection criteria and fake 
an electron. From the shape of the spectrum for these 
events we determine how likely it is for these events to 
have sufficient to enter our W sample. 

We determine this shape by selecting isolated electro- 
magnetic clusters that have > 200 and (Ttrk > 10. 
Almost all electrons fail this cut, so that the remaining 



sample consists almost entirely of hadrons. We use data 
taken by a trigger without the j^t requirement to study 
the efficiencies of this cut for jets. For j^t < 10 GeV 
we find 1973 such events, while in the same sample 3674 
satisfy our electron selection criteria. If we normalize the 
background spectrum to the electron sample we obtain 
an estimate of the hadronic background in an electron 
candidate sample. Figure |5j shows the spectra of 
both samples, normalized for j^t < 10 GeV. 

In the data collected with the W trigger we find 204 
events that satisfy all the fiducial and kinematic cuts. 



listed in Sec. [V for the W sample, and have x > 200 



and CTtrk > 10. We therefore estimate that 374 back- 
ground events entered the signal sample. This corre- 
sponds to a fraction of the total W sample after all cuts 
of /had = (1.3 ± 0.2) %. For a looser cut on the recoil px, 
UT < 30 GeV, we find /had = (1-6 ± 0.3) %. The error is 
dominated by uncertainty in the relative normalization 
of the two samples at low ^t- Figure |5^ shows the back- 
ground fraction as a function of luminosity. There is no 
evidence for a significant luminosity dependence. We use 
the background events with pt{v) > 25 GeV to estimate 
the shape of the background contributions to the pT(e), 
Pt{i'), and my spectra (Fig. [5q). 



C. Z ^ee 

To estimate the fraction of Z ^ ee events which sat- 
isfy the W selection, we use a Monte Carlo sample of 
approximately 100,000 Z — > ee events generated with 
the HERWIG program and a detector simulation based on 
GEANT. The boson px spectrum generated by herwig 
agrees reasonably well with the calculation in Ref. Q. 
Z — > ee decays typically enter the W sample when one 
electron satisfies the W cuts and the second electron is 
lost or mismeasured, causing the event to have large j^t- 

Approximately 1.1% of the Z ^ ee events have an 
electron with pseudorapidity > 4.0, which is the 
acceptance limit of the end calorimeters. The frac- 
tion of Z — > ee events which contain one electron with 
|77(ei)| < 1.0 and pT(e) > 25 GeV, and another with 
177(62)1 > 4.0 is approximately 0.04%. The contribution 
from the case of an electron lost through the beampipe 
is therefore relatively small. 

An electron is most frequently mismeasured when it 
goes into the regions between the CC and one of the 
ECs, which are covered only by the hadronic section of 
the calorimeter. These electrons therefore can not be 
identified and their energy is measured in the hadronic 
calorimeter. Large is more likely for these events 
than when both electrons hit the EM calorimeters. The 
mismeasured electron contributes to the recoil when the 
event is treated as a M^. The fraction of Z events in the 
W sample therefore depends on the ut cut. 

We find that 10,987 Monte Carlo events pass the CC- 
CC Z ^ ee selection, and 758 (1,318) pass the W selec- 
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tion with a recoil cut of 15 (30) GeV. The fraction of Z 
events in the W sample is therefore jz — (0.42 ± 0.08) % 
for UT < 15 GeV and (0.62 ± 0.08) % for ut < 30 GeV. 
The uncertainties quoted include systematic uncertain- 
ties in the matching of momentum scales between Monte 
Carlo and collider data. Figure |5^ shows the distribu- 
tions of pt{&), Pt{v), and tot for the events that satisfy 
the W selection. 



data spectra we add the background contributions with 
the shapes and normalizations described in Sec. IX to the 
signal spectra. We normalize the spectra within the fit 
interval and interpret them as probability density func- 
tions to compute the likelihood 



N 



L{m) = JJpi(TO)'' 



(51) 



D. W 



hadrons + X 



We estimate the background due to W ^ tv followed 
by a hadronic tau decay based on two Monte Carlo sam- 
ples. In a sample oiW ~* tv hadrons + X simulated 
using GEANT, 65 out of 4,514 events pass the fiducial 
and kinematic cuts of the W sample. We use a sam- 
ple of W ^ TV ^ hadrons -f X simulated by replacing 
the electron shower \nW ^ ev decays from collider data 
with the hadrons from a tau decay, generated by a Monte 
Carlo simulation, to estimate the probability of the tau 
decay products to fake an electron. Of 552 events that 
pass the fiducial and kinematic cuts 145 pass the elec- 
tron identification criteria. With the hadronic branching 
fraction for taus, J5(r — > hadrons) = 64% we estimate a 
contamination of the W sample of 0.24% from hadronic 
tau decays. The expected background shapes are plotted 
in Fig. 



E. Cosmic Rays 

Cosmic ray muons can cause backgrounds when they 
coincide with a beam crossing and radiate a photon of 
sufficient energy to mimic the signature of the electron 
from W ^ ev decays. We measure this background by 
searching for muons near the electrons in the W signal 
sample. The muons have to be within 10° of the electron 
in azimuth. Using muon selection criteria similar to those 
in Ref. |Q we observe 18 events with such muons in the 
W sample. We estimate the fraction of cosmic ray events 
in the W sample to be 0.2 ± 0.1%. The effect of this 
background on the W mass measurement is negligible. 



X. MASS FITS 

A. Maximum Likelihood Fitting Procedure 

We use a binned maximum likelihood fit to extract the 
W mass. Using the fast Monte Carlo program we com- 
pute the rriT, pxie), and pt{v) spectra for 200 hypothe- 
sized values of the W mass between 79.4 and 81.4 GeV. 
For the tot spectrum we use 100 MeV bins and for the 
lepton pt spectra we use 50 MeV bins. The statistical 
precision of the spectra for the W mass fit corresponds 
to about 4 million W decays. When fitting the collider 



where Pi{m) is the probability density for bin i, assuming 
Mw = TO, and is the number of data entries in bin i. 
The product runs over all N bins inside the fit interval. 
We fit — ln(L(TO)) with a quadratic function of m. The 
value of m at which the function assumes its minimum 
is the fitted value of the W mass and the 68% confidence 
level interval is the interval in m for which — ln(L(TO)) is 
within half a unit of its minimum. 

As a consistency check of the fitting procedure we gen- 
erate 105 Monte Carlo ensembles of 28,323 events each 
with Mvi/=80.4 GeV. We then fit these ensembles with 
the same probability density functions as the collider 
data spectra, except that we do not include the back- 
ground contributions. Table VI lists the mean, rms, and 



correlation matrix of the fitted values. 



B. Electron px Spectrum 

We fit the prie) spectrum in the region 30 < prie) < 
50 GeV. There are 22,898 events in this interval. The 
data points in Fig. ^ represent the prie) spectrum from 
the W sample. The solid line shows the sum of the simu- 
lated W signal and the estimated background for the best 
fit, and the shaded region indicates the sum of the esti- 
mated hadronic, Z ee, and W tv hadrons + X 
backgrounds. The maximum likelihood fit gives 



Mw = 80.475 ± 0.087 GeV 



(52) 



for the W mass. 

As a goodness-of-fit test, we divide the fit interval into 
0.5 GeV bins, normalize the integral of the probability 
density function to the number of events in the fit inter- 
val, and compute = J2iLi iVi ~ PiY IVi- The sum runs 
over all N bins, yi is the observed number of events in 
bin i, and Pi is the integral of the normalized probability 
density function over bin i. The parent distribution is 
the distribution for TV — 2 degrees of freedom. For 
the spectra in Fig. |5^ we compute — 40.6. For 40 
bins there is a 35% probability for ^ 40.6. Figure |58| 
shows the contributions Xi = [Ui ^ Pi)/ y/ui to for the 
40 bins in the fit interval. 

We also compare the observed spectrum to the prob- 
ability density function using the Kolmogorov-Smirnov 
test. For a comparison within the fit window we obtain 
K = 0.81 and for the entire histogram k = 0.83. 

Figure ^ shows the sensitivity of the fitted mass value 
to the choice of fit interval. The points in the two plots 
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indicate the observed deviation of the fitted mass from 
the value given in Eq. We expect some variation due 
to statistical fluctuations in the spectrum and systematic 
uncertainties in the probability density functions. We 
estimate the effect due to statistical fluctuations using 
the Monte Carlo ensembles described above. We expect 
the fitted values to be inside the shaded regions indicated 
in the two plots with 68% probability. The dashed lines 
indicate the statistical error for the nominal fit. 

All tests show that the probability density function 
provides a good description of the observed spectrum. 



within the fit window and k, = 0.69 for the entire his- 
togram. Figure ^ shows the deviation x between data 
and fit. There is an indication of a systematic deviation 
between the observed spectrum and the resolution func- 
tion. This effect is not very significant. For example, 
when we increase the hadronic resolution parameter amb 
in the simulation to 1.11, which corresponds to about 
1.5 standard deviations, this indication of a deviation 
between data and Monte Carlo vanishes. 



B. Luminosity Dependence 



C. Transverse Mass Spectrum 

Figure ^ shows the rriT spectrum. The points are the 
observed spectrum, the solid line shows signal plus back- 
ground for the best fit, and the shaded region indicates 
the estimated background contamination. We fit in the 
interval 60 < rriT < 90 GeV. There are 23,068 events in 
this interval. Figure |6l| shows —\n{L{m)/Lo) for this fit 
where Lq is an arbitrary number. The best fit occurs for 



We divide the W and Z data samples into four lumi- 



nosity bins 



Mw = 80.438 ± 0.070 GeV. 



(53) 



Figure ^ shows the deviation of the data from the fit. 
Summing over all bins in the fitting window, we get = 
79.5 for 60 bins. For 60 bins there is a 3% probability 
to obtain a larger value. The Kolmogorov-Smirnov test 
gives K = 0.25 within the fit window and k = 0.84 for 
the entire histogram. Figure ^ shows the sensitivity of 
the fitted mass to the choice of fit interval. 

In spite of the somewhat large value of there is no 
structure apparent in Fig. ^ that would indicate that 
there is a systematic difference between the shapes of the 
observed spectrum and the probability density function. 
The large can be attributed to a few bins that are 
scattered over the entire fit interval, indicating statistical 
fluctuations in the data. This is consistent with the good 
Kolmogorov-Smirnov probability which is more sensitive 
to the shape of the distribution and insensitive to the 
binning. 



XI. CONSISTENCY CHECKS 

A. Neutrino pr Spectrum 

As a consistency check, we also fit the pt{v) spectrum, 
although this measurement is subject to much larger sys- 
tematic uncertainties than the tjit and pt{g) fits. Fig- 
ure ^ shows the observed spectrum (points) , signal plus 
background for the best fit (solid line), and the estimated 
background (shaded region). For the fit interval 30 < 
Pt{v) < 50 GeV the fitted mass is Mw = 80.37 ± 0.11 
GeV, in good agreement with the niT and pxie) fits. We 
compute = 31.8. The probability for a larger value 
is 75%. The Kolmogorov-Smirnov test gives k = 0.20 



£ < 5 X lO^Ocm-^s-i, 



30^ 



5 X lO'^'J < £ < 7 X 10 
7 X 10^° < £ < 9 X 103°cm-2s-i, 



and generate resolution functions for the luminosity dis- 
tribution of these four subsamples. We fit the transverse 
mass and lepton px spectra from the W samples and the 
dielectron invariant mass spectra from the Z samples in 
each bin. The fitted masses are plotted in Fig. The 
errors are statistical only. We compute the x^ with re- 
spect to the W mass fit to the mr spectrum from the 
entire data sample. The x^ per degree of freedom (dof) 
for the pT{e) fit is 1.9/4 and for the pt{i^) fit is 2.4/4. 
The rriT fit has a x^/dof of 2.7/3. The solid and dashed 
lines in the top plot indicate the W mass value and sta- 
tistical uncertainty from the fit to the ttit spectrum of 
the entire data sample. All measurements are in very 
good agreement with this value. In the bottom plot the 
lines indicate the Z mass fit to the m(ee) spectrum of 
the entire Z data sample. The measurements in the four 
luminosity bins have a x^/dof of 1.0/3. 



C. Dependence on ut Cut 

We change the cuts on the recoil momentum ut and 
study how well the fast Monte Carlo simulation repro- 
duces the variations in the spectra. We split the W sam- 
ple in two subsamples with M|| > and U|| < 0. In the 
simulation we fix the W mass to the value from the m-p 
fit in Eq. Figures |6^^ show the itlt, prie), and 
Pt{i^) spectra from the collider data for the subsamples 
with u\\ > and U|| < and the corresponding Monte 
Carlo predictions. Table VII lists the results of com- 



parisons of collider data and Monte Carlo spectra using 
the Kolmogorov-Smirnov test. Although there is signif- 
icant variation among the shapes of the spectra for the 
different cuts, the fast Monte Carlo models them well. 



Table VII also lists the results of comparisons of collider 
data and Monte Carlo spectra for a W sample selected 
with Ut < 30 GeV which consists of 32,361 events. 



20 



D. Dependence on Fiducial Cuts 

We divide the azimuth of the recoil momentum, 4'{R), 
into eight bins. This binning is sensitive to azimuthal 
nonuniformities in the recoil momentum measurement, 
e.g. because of background from the Main Ring. Fig- 
ure ^ shows the fitted W mass values versus 4>{R) ■ The 
Main Ring is located at ~ 7r/2 and any biases caused by 
background from the Main Ring should appear as struc- 
ture in this direction or in the opposite direction. The 
rms of the eight data points is 124 MeV, consistent with 
the statistical uncertainty of 200 MeV for the data points. 
Thus the data are consistent with azimuthal uniformity. 

We divide the azimuthal direction of the electron, 0(e), 
into 32 bins corresponding to the 32 azimuthal modules 
of the CC-EM. Figure ^ shows the fitted W mass val- 
ues versus 4>{e). The statistical uncertainty of the data 
points is 400 MeV and the rms of the 32 points is 600 
MeV. Thus there is a 0.6% nonuniformity in the response 
of the CC-EM, consistent with the module-to-module cal- 
ibration of 0.5% ||. 

Finally, we fit the rriT spectrum from the W sample 
and the m(ee) spectrum from the Z sample for different 
pseudorapidity cuts on the electron direction. We use 
cuts of |7/(e)| < 1.0, 0.7, 0.5, and 0.3. Figure ^ shows 
the change in the W mass versus the ri{e) cut using the 
electron energy scale calibration from the corresponding 
Z sample. The shaded region indicates the statistical 
error. Within the uncertainties the mass is independent 
of the 77(e) cut. 



XII. UNCERTAINTIES IN THE 
MEASUREMENT 

A. Statistical Uncertainties 



Table VIII lists the uncertainties, rounded to the near- 



est 5 MeV, in the W measurement due to the finite sizes 
of the W and Z samples used in the fits to the my, Pt{s), 
Pt{v), and m{ee) spectra. The statistical uncertainty 
due to the finite Z sample propagates into the W mass 
measurement through the electron energy scale aEM- 



B. W Production and Decay Model 

1. Sources of Uncertainty 

Uncertainties in the W production and decay model 
arise from the following sources: the phenomcnological 
parameters in the calculation of the pt{W) spectrum, the 
choice of parton distribution functions, radiative decays, 
and the W boson width. In the following we describe 
how we assess the size of the systematic uncertainties 
introduced by each of these. We summarize the size of 



the uncertainties in Table IX, rounded to the nearest 5 
MeV. 



2. W Boson Pt Spectrum 



In Sec. VIII we determine 32 so that the predicted 
priee) spectrum agrees with the Z data. In order to 
quantify the uncertainty in the boson pT spectra, we need 
to consider variations in all four parameters, Aqcd, 5i, 
52, and (73. We use a series of modified CTEQ3M parton 
distribution functions fit with Aqcd fixed at discrete val- 
ues to study the variations in the pT{ee) spectrum 
and the fitted W boson mass with these parameters. 

We cannot constrain all these four parameters simul- 
taneously by using only our Z data. We therefore in- 
troduce an external constraint on Aqcd- The CTEQ3M 
fits prefer Aqcd = 158 MeV but are also consistent with 
somewhat higher values . Other measurements give a 
combined value of Aqcd = '209^11 MeV ^ . All data are 
consistent with Aqcd between 150 and 250 MeV, which 
we use as the range over which Aqcd is allowed to vary. 

The requirement that the fast Monte Carlo prediction 
for the average priee) over the range pT(ee) < 15 GeV, 
corrected for background contributions, must agree with 
the value observed in the Z data, {priee)) = 6.05 ± 0.07 
GeV, couples the values of Aqcd and 32- Figure [T^ shows 
a plot of g2 versus Aqcd ■ For any pair of values on the 
curve the fast Monte Carlo predicts a value of {priee)) 
that agrees with the Z data. For any fixed value of Aqcd , 
52 is determined to a precision of 0.12 GeV'^. This error 
includes the statistical uncertainty (0.09 GeV^) and the 
systematic uncertainty due to normalization and shape 
of the background (0.07 GeV^). All other uncertainties, 
e.g. due to electron momentum resolution and response 
or selection biases, are negligible. 

If we fix Aqcd and 32, the requirement that the aver- 
age Priee) predicted by the fast Monte Carlo agree with 
the data allows an additional variation in the parameters 
gi and g^. The residual uncertainty in the measured W 
boson mass due to this variation, however, is small com- 
pared to the uncertainty due to the variation allowed in 
52 and Aqcd and we neglect it. Finally, we obtain the 
uncertainties in the fitted W boson mass listed in Ta- 



ble IX 



3. Parton Distribution Functions 

The choice of parton distribution function used to de- 
scribe the momentum distribution of the constituents of 
the proton and antiproton affects several components of 
the model: the parton luminosity slope /?, and the rapid- 
ity and transverse momentum spectrum of the W. 

Using several modern parton distribution function sets 
as input to the fast Monte Carlo model, we generate niT 
and lepton px spectra. In each case we use the value of g2 
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measured for that part on distribution function set using 
our Z data (Sec. VIII ). We then fit them in the same 
way as the spectra from coUider data, i.e. using MRS A' 
parton distribution functions. Table |X| Usts the variation 
of the fitted W mass values relative to MRS A'. 

The MRSA' and CTEQ3M parton distribution func- 
tions use the measured W charge asymmetry in pp colli- 
sions as input to the fit. MRSD-' and CTEQ2M do 
not explicitly use the asymmetry. The asymmetry pre- 
dicted by MRSD— ' agrees with the measurement; that 
of CTEQ2M disagrees at the level of four standard de- 
viations. We include CTEQ2M in our estimate of the 
uncertainty to provide an estimate of the possible vari- 
ations with a rather large deviation from the measured 
asymmetry. 



J^. Parton Luminosity 



The uncertainty of 10^'^ GeV^^ in the parton luminos- 
ity slope /3 (Sec. [V|) translates into an uncertainty in the 
fitted W mass. We estimate the sensitivity in the fitted 
W mass by fitting Monte Carlo spectra generated with 
different values of /3. 



used to fit the mass in the final state. Finally, the cal- 
culation includes only processes in which a single photon 
is radiated. We use the code provided by the authors 
of Ref. [|l] to estimate the shift introduced in the mea- 
sured W mass by neglecting two-photon emission. We 
find that two photons, with -pT > 100 MeV and sepa- 
rated by AR > 0.3 from the electron, are radiated in 
about 0.24% of all W ^ ev decays. This reduces the 
mean value of my within the fit window by 3 MeV. In 
1.1% of all Z ^ ee decays two photons, with pt > 100 
MeV and separated by AR > 0.3 from the electrons, are 
radiated. We add the dielectron mass spectrum of these 
Z 6677 events to our simulated Z boson lineshape and 
fit the modified lineshape. The fitted mass decreases by 
10 MeV. This shift requires an adjustment of the energy 
scale calibration factor a^M by 10^**. Neglecting two- 
photon emission in both W and Z boson decays then 
increases the measured W mass by about 5 MeV. Since 
this effect is an order of magnitude smaller than the sta- 
tistical uncertainty in our measurement we do not correct 
for it, but add it in quadrature to the uncertainty due to 
radiative corrections. 



6. W Boson Width 



5. Radiative Decays 

We assign an error to the modeling of radiative de- 
cays based on varying the detector parameters Eq and 
i?o (Sec. ^). Eq defines the minimum photon energy 
generated and corresponds to a cut-off below which the 
photon does not reach the calorimeter. Rq defines the 
maximum separation between the photon and electron 
directions above which the photon energy is not included 
in the electron shower. In general, radiation shifts the 
fitted mass down for the transverse mass and electron 
fits, because for a fraction of the events the photon en- 
ergy is subtracted from the electron. Hence increasing 
Rq decreases the radiative shift. Similarly, decreasing Eq 
decreases the radiative shift. Both the fitted W and Z 



masses depend on these parameters. Table XI lists the 



change in the fitted masses if radiative effects are turned 
off completely. To estimate the systematic error, we fit 
Monte Carlo spectra generated with different values for 
Eq and Rq . For the low value of Eq = 50 MeV that we use 
in the simulation, the dependence of the fits on this pa- 
rameter is negligible. The changes in the mass fits when 
varying Rq by ±0.1 are also listed in Table XI. After 



propagating the change in the Z mass into the electron 
response the result of the W mass measurement changes 
by about 15 MeV for all three spectra. 

There are also theoretical uncertainties in the radia- 
tive decay calculation. Initial state QED radiation is not 
included in the calculation of Ref. However, initial 
state radiation does not affect the kinematic distributions 



To determine the sensitivity of the fitted W mass to the 
W width, we generate mr and lepton px spectra using 
the fast Monte Carlo model with a range of widths and 
fit them with the nominal templates. The uncertainty on 
the fitted W mass correspond to the uncertainty in the 
measured value oiTw = 2.062±0.059 GeV M. 



C. Detector Model Parameters 

The uncertainties on the p aram eters of the detector 
model determined in Sees. VI -VII translate into uncer- 



tainties in the W mass measurement. We study the sen- 
sitivity of the W mass measurement to the values of the 
parameters by fitting the data with spectra generated by 
the fast Monte Carlo with modified input parameters. 

Table XII lists the uncertainties in the measured W 
mass, caused by the individual parameters. We assign 
sets of correlated parameters to the same item in the ta- 
ble. Correlations between items are negligible. For each 
item the uncertainty is determined to typically 5 MeV 
for the niT fit and 10 MeV for the lepton px fits. We 
therefore round them to the nearest 5 MeV in the table. 
To achieve this precision 10-20 million W ev decays 
are simulated for each item. 

The residual calorimeter nonlinearity is parametrized 
by the offset Jem- Calorimeter uniformity refers to a pos- 
sible nonuniformity in response as a function of rj. It is 
limited by the test beam data |l^ . The electron momen- 
tum resolution is parametrized by cem- The electron 
angle calibration includes the effects of the parameters 
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acDC and acc, discussed in Appendices ^ and ^ The 
recoil resolution is parametrized by Q;,„b and Sicc and the 
response by arcc and Prcc- Electron removal refers to the 
bias Au|| introduced in the U|| measurement by the re- 
moval of the cells occupied by the electron. Selection bias 
refers to the U|| efficiency. 



D. Backgrounds 

We determine the sensitivity of the fit results to the as- 
sumed background normalizations and shapes by repeat- 
ing the fits to the data wit h var ied background shapes 
and normalizations. Tabic XllI lists the uncertainties 
rounded to the nearest 5 MeV. 



XIII. RESULTS 

We present a precision measurement of the mass of the 
W boson. From a fit to the transverse mass spectrum, 
we measure 



Mw = 80.44 ± O.lO(stat) ± 0.07{syst) GeV. 



(54) 



Adding all errors in quadrature gives 115 MeV. Since we 
calibrate the electron energy scale against the known Z 
mass, we effectively measure the W and Z mass ratio 



w 



Mz 



0.8821 ± O.OOll(stat) ± 0.0008(s?/st). 



(55) 



A fit to the transverse momentum spectrum of the decay 
electrons gives 



M, 



w 



80.48 ± O.ll(siat) ± Q.m{syst) GeV. 



(56) 



Adding all errors in quadrature gives 140 MeV. As ex- 
pected, the measurement from the mr spectrum has a 
larger uncertainty from detector effects (65 MeV) than 
that from the prie) spectrum (50 MeV). On the other 
hand the tot fit is less sensitive to the W production 
model (30 McV) than the prie) fit (75 MeV). The good 
agreement between the two results indicates that we un- 
derstand the ingredients of our model and their uncer- 
tainties. In the end, the to^ fit gives the more precise 
result and we quote this as our final result. However the 
fit to the pt (e) spectrum may become more competitive 
in the future with larger data samples and better con- 
straints on th e W production dynamics. 

Table XIV lists the D0 W mass measurements from 
fits to the Tot spectra from the 1992-1993 |0 and the 
1994-19 95 dat a sets and their uncertainties. As indicated 
in Table XIV, some errors are common to the two mea- 



surements. Since both analyses use the same W produc- 
tion and decay model we assign the uncertainties quoted 



in Sec. XII B| to both measurements. The precision of 
the electron angle calibration has improved compared to 
Rcf. and we use the reduced uncertainty for both 
measurements. All uncertainties due to detector model 
parameters, which were measured using statistically in- 
dependent data sets, are uncorrelated because their pre- 
cision is dominated by statistical fluctuations. In order 
to combine the two measurements we weight them by 
their uncorrelated errors Sa and 



M, 



w 



Ma/51 



V5l 



The uncertainty is then given by 



5Mw = 



^1 51 + 1/51 



5\ 



(57) 



(58) 



where 5 i s the common uncertainty from the third column 
in Table XIV. The combination of the D0 measurements 



from the 1992-1993 and 1994-1995 data gives 



w 



80.43 ±0.11 GeV. 



(59) 



The D0 measurement is in good agreement with pre- 
vious measurements and is more precise than all the pre- 
viously published measurements combined. Table XV 



lists previously published measurements with uncertain- 
ties below 500 MeV. A global fit to all electroweak mea- 
surements from the LEP experiments predicts Mw = 
80.278 ± 0.049 GeV §. Figure ^ gives a graphical rep- 
resentation of these data. 

We evaluate the radiative corrections Ar, defined in 
Eq. |l|. Our measurement of Adw from Eq. ^ leads to 



Ar = -0.0288 ±0.0070, 



(60) 



In 



4.1 standard deviations from the tree level value 
Fig. 75 we compare the measured W and top quark 
the values predicted by the Standard 

Also shown 
l2ltl for a 



masses |gO|l to 

Model for a range of Higgs mass values 
is the prediction from the calculation in Ref. 
model involving supersymmetric particles assuming the 
chargino, Higgs, and left-handed selectron masses are 
greater than 90 GeV. The measured values are in agree- 
ment with the prediction of the Standard Model. 
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that a scale of acDC = 0.989 ±0.001 is needed to correct 
the track. 

Combining all measurements of acDC gives acDC = 
0.988 ± 0.001, which we use in the reconstruction of the 
electrons in the W and Z data samples. 



APPENDIX A: TRACK POSITION 
CALIBRATION 

We use cosmic ray muons which traverse the entire de- 
tector and pass close to the beam position to calibrate 
the z-measurement of the track in the CDC. We predict 
the trajectory of the muon through the central detec- 
tor by connecting the incoming and outgoing hits in the 
innermost muon chambers by a straight line. The cen- 
ter of gravity of the incoming and outgoing CDC tracks 
are then calibrated relative to this line. Figure ^ shows 
the difference between the predicted and the actual z- 
positions of the track centers of gravity. These data are 
fit to a straight line. We find the track position must be 
scaled by the fitted slope, acDC = 0.9868 ± 0.0004. 

We also use a sample of low-pT dimuon events from 
pp collisions where both muons originate from the same 
interaction vertex. We reconstruct the muon trajectories 
from their hits in the innermost muon chambers and the 
CDC. For both muons we determine the point of closest 
approach of the trajectory to the beam, Zvtx(A*)- We then 
scale the z-position of the CDC track to minimize 



events 



^vtx ^vtx 



(Al) 



where cr^ is chosen so that the minimum value of 
equals the number of events minus one. The minimum 
occurs at acoc = 0.9863±0.0011. The same analysis ap- 
plied to a Z — > /i^ sample gives acDC = 0.9878 ± 0.0014 
and is shown in Fig. ^ 

A scintillating fiber detector was inserted between the 
CDC and the CC to calibrate the track z-position. The 
detector is built from 20 modules, each constructed on 
an aluminum support plate 93.4 cm long and 16.5 cm 
wide. Scintillating fibers, 12.7 cm long, were laid across 
the width of the module every 11.43 cm along the support 
plate. The eight scintillating fibers on each module were 
connected to a clear waveguide and read out with a pho- 
tomultiplier tube. The modules are mounted lengthwise 
along the cylinder of the CDC with half of the modules 
covering +z and the other half —z. In the r-(j) view each 
module subtends 7r/16 radians with the fibers running az- 
imuthally. Because of spacial constraints not the entire 
CDC was covered. 

When a fiber is hit by a charged particle the z-position 
of the associated track, at the fiber radius, is compared 
with the fiber z-position. The z-position of the track at 
the radial position of the fiber is determined from the di- 
rection and center of gravity of the track. By comparing 
the z-position of the track and the hit fiber, we determine 



APPENDIX B: 



ELECTRON SHOWER POSITION 
ALGORITHM 



We determine the position of the electron shower cen- 
troid Xcai — {xced, ycai, ^cai) in the calorimeter from the 
energy depositions in the third EM layer by computing 
the weighted mean of the positions Xi of the cell centers, 



WiXi 



Wi 



The weights are given by 



Wi 



max 0, 1^0 + log 



E{e) 



(Bl) 



(B2) 



where Ei is the energy in cell «, is a parameter which 
depends upon ri{e), and E(e) is the energy of the electron. 
We calibrate the algorithm using Monte Carlo electrons 
simulated using GEANT and electrons from the Z ^ ee 
data. We apply a polynomial correction as a function of 
Zeal and 9{e) based on the Monte Carlo electrons. We 
refine the calibration with the Z ^ ee data by exploit- 
ing the fact that both electrons originate from the same 
vertex. Using the algorithm given by Eq. Al we deter- 



mine a vertex for each electron from the shower centroid 
and the track center of gravity. We minimize the dif- 
ference between the two vertex positions as a function 
of a scale factor acc- More complex correction func- 
tions do not improve the . The correction factor is 
acc — 0.9980 ±0.0005, where the error includes possible 
variations of the functional form of the correction. 
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FIG. 1. Loop diagrams contributing to the W boson mass. 
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FIG. 3. The mr spectrum for W bosons with = 

( ), with the correct qr distribution (•), and with detector 

resolutions (shaded). 
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FIG. 5. A cutaway view of the D0 calorimeter and track- 
ing system. 
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FIG. 6. The pedestal spectrum of a central calorime- 
ter cell, where the mean pedestal has been subtracted. The 
shaded region are the events removed by the zero-suppression. 
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FIG. 7. The response of the liquid argon in the central 
calorimeter as monitored by a and /3 sources. 
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FIG. 8. The percentage change in the central calorimeter 
stains over the course of the run. 
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FIG. 9. The change in the central calorimeter pedestals 
over the course of the run. 
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FIG. 10. The fractional deviation of the reconstructed 
electron energy from the beam momentum from beam tests 
of a CC-EM module. 
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FIG. 12. The relative eflSciency of the Level 2 electron 
filter for a threshold of 20 GcV. The arrow indicates the cut 
applied in the final event selection. 




FIG. 11. The fractional electron energy resolution mea- 
sured in beam tests of a CC-EM module for the data (•) and 

the parameterization ( ). 



FIG. 13. The efficiency of a 15 GeV Level 2 require- 
ment. The arrow indicates the cut applied in the final event 
selection. 
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FIG. 14. Illustration of momentum vectors in the trans- 
verse plane for Z ^ ee candidates. The vectors drawn with piQ ig. Distributions of the electron identification vari- 
thick Unes are directly measured. ables. The arrows indicate the cut values. 
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FIG. 18. Polarization of the W produced in pp collisions 
if the quark comes from the proton (left) and if the antiquark 
comes from the proton (right). The thick arrows indicate the 
orientation of the particle spins. 
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FIG. 19. The calculations of ai( ) and 02 ( ) as a 

function of the transverse momentum of the W boson. 
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FIG. 20. The distribution of AR{e'y) of photons from 
W — > evy decays that are reconstructed as separate objects 
(shaded) and those that are not, either because they are too 
close to the electron or too low in energy ( ). 




FIG. 21. The distribution of 2vtx(ei) — Zvtx(e2) for 
the Z —f ee sample (•) and the fast Monte Carlo simula- 
tion ( ). 
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FIG. 22. The transverse energy flow into 1x5 tower seg- 
ments as a function of azimuthal separation from the electron 
in the W sample. 
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. 24. The luminosity dependence of (A'U||). 
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FIG. 25. The variation of (Am||) as a function of M||. The 
region between the arrows is populated by the W sample. 
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FIG. 23. The distribution of Awm in the W signal sample. 
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FIG. 26. The transverse energy flow into 1x5 tower seg- 
ments as a function of the azimuthal separation from the elec- 
tron for the electrons from W ^ ev decays (•) and the su- 
perimposed Monte Carlo electron sample ( ). 



31 




FIG. 27. The isolation spectrum for five different re- 
gions, it|| < -15, -15 < U|| < -5, -5 < U|| < 5, 5 < ity < 15, pjQ 29. The electron selection efhciency as a function of 

> 15 GeV (from top to bottom), for the electrons from ^ 
ly — > ez/ decays (•) and the superimposed electron sample 
(shaded) . 
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FIG. 28. The mean isolation versus My for the W elec- 
tron sample (o) and the superimposed Monte Carlo electron 
sample (•). 
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FIG. 30. The background-subtracted distribution. 
The superimposed curve shows the Monte Carlo simulation. 
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FIG. 31. The dielectron invariant mass spectrum for the 
J/'tp ee sample (histogram) and the background sample 
(•). The smooth curve is a fit to the data. 
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FIG. 33. The distribution of m(ee) versus fz for the 
CC-CC Z ~*ee sample. 
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FIG. 32. 
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The dielectron mass spectrum for the back- 
sample to the CC/CC Z sample. The fit is an 
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FIG. 34. The 68% confidence level contours in qem and 
Sem from the J/iIj, tv", and Z data. The inset shows an ex- 
panded view of the region where the is minimized. 
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FIG. 35. The dielectron mass spectrum from the CC-CC 
Z sample. The superimposed curve shows the maximum like- 
lihood fit and the shaded region the fitted background. 



34 




70 80 90 100 110 
m(ee) (GeV) 



70 80 90 100 110 
m(ee) (GeV) 



FIG. 36. The dielectron mass spectra from (a) the 
tight/tight and (b) the tight/loose CC-CC Z samples. The 
curves show the fitted Monte Carlo spectra. 
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FIG. 37. The dielectron mass resolution versus the con- 
stant term cem- 
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FIG. 38. The recoil momentum response in the Monte 
Carlo Z sample as a function of qr- 
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FIG. 39. The average Prj{ee) + versus Pr)(ee) for the Z 
data (•) and the fast Monte Carlo simulation (o) . 
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FIG. 41. The width of the ?7-balance distribution versus 
Pr, (ee) for the Z data (•) and the fast Monte Carlo simulation 
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FIG. 43. The r;-balance distribution for the Z data (•) 
and the fast Monte Carlo simulation ( ). 
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FIG. 44. The ^-balance distribution for the Z data (• 
and the fast Monte Carlo simulation ( ). 
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FIG. 46. The resolution for transverse momentum bal- 
ance, (tt, versus the transverse energy flow, St, for minimum 
bias events (•). The smooth curve is a fit (Eq. Eh. 
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FIG. 45. The transverse energy flow in the W (•) and Z 
( ) data. 
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FIG. 47. The spectrum for the W data (•) and the 
Monte Carlo simulation ( ). 
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FIG. 48. The u± spectrum for the W data (•) and the 
Monte Carlo simulation ( ). 
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FIG. 50. The azimuthal difFcrcnco between electron and 
recoil directions for the W data (•) and the Monte Carlo 
simulation ( ). 
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FIG. 49. The recoil momentum (ut) spectrum for the W 

data (•) and the Monte Carlo simulation ( ). The arrow 

shows the location of the cut. 
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FIG. 51. Comparison of the pT(ee) data (•) and simula- 
tion ( ) for the best fit g2 using MRSA' parton distribution 

functions. 



38 




FIG. 52. The background parameterizations for the 
priee) spectrum. 
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FIG. 53. The predicted priee) spectra after detector 
simulation using MRSA' parton distribution functions and 
52 = 0.18, 0.58, and 0.98 GeV^ 
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FIG. 55. The fraction of hadron background as a function 
of luminosity. 
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FIG. 54. The spectra of a sample of events passing elec- 
tron identification cuts ( ) and a sample of events failing 

the cuts (•). 
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FIG. 56. Shape of rriT, pri^e), and priy) spectra from 

hadron ( ), Z ( ), and t ^hadron backgrounds ( ) 

with the proper relative normalization. 
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FIG. 57. Spectrum of pT{e) from the W data. The su- 
perimposed curve shows the maximum likelihood fit and the 
shaded region the estimated background. 
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FIG. 59. Variation of the fitted mass with the pT(e) fit 
window limits. See text for details. 
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FIG. 58. The x distribution for the fit to the prie) spec- 
trum. 
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FIG. 60. Spectrum of rriT from the W data. The superim- 
posed curve shows the maximum likelihood fit and the shaded 
region shows the estimated background. 
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FIG. 61. The likelihood function for the mr fit. 



80 82.5 85 87.5 90 92.5 95 97.5 100 

upper window limit (GeV) 

FIG. 63. Variation of the fitted mass with the mx fit 
window limits. See text for details. 
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The X distribution for the fit to the rriT spec- 
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FIG. 64. Spectrum of pt{v) from the W data. The su- 
perimposed curve shows the maximum likelihood fit and the 
shaded region shows the estimated background. 
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FIG. 65. The x distribution for the fit to the Pt(i') spec- 
trum. 
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FIG. 67. Spectra of mr from W data with u\\ < (o) and 
W|| > (•) compared to Monte Carlo simulations ( ). 
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FIG. 66. The fitted W boson masses (a) in bins of lu- 
minosity from the my (•), Prie) (o), and prii^) (*) fits (the 
points are ofiiset for clarity) and the fitted Z boson masses 
(b). The solid line is the central value for the mr and m(ee) 
mass fits respectively over the entire luminosity range and the 
dashed lines are the statistical errors. 
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FIG. 68. Spectra of prie) from W data with wy < (o) 
and W|| > (•) compared to Monte Carlo simulations ( ). 
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FIG. 69. Spectra of Pt{i^) from W data with U|| < 
and W|i > (•) compared to Monte Carlo simulations ( — 
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FIG. 70. The variation in the W mass from the nir fit as 0.2 0.4 0.6 0.8 1 

a function of <j>{R). r|(e) Cut 



FIG. 72. The variation in the W mass versus the 77(e) cut. 
The shaded region is the expected statistical variation. 




FIG. 73. Value of g2 as a function of Aqcd- The error bar 
indicates the uncertainty in g2 for fixed Aqcd . 
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FIG. 74. A comparison of this measurement with previ- 
ously published W mass measurements (Table XV). The 
shaded region indicates the predicted W mass value from 
global fits to the Z lineshape data j^. 
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FIG. 75. A comparison of the W and top quark mass mea- 
surements by the D0 collaboration with the Standard Model 
predictions for different Higgs boson masses [^. The width 
of the bands for each Higgs boson mass value indicates the 
uncertainty due to the error in a(Mf ). Also shown is the 
range allowed by the MSSM H]. 
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FIG. 77. The x versus qcdc value. The arrows indicate 
the statistical error on the fit. 

TABLE I. Electron selection criteria. 

variable CC (loose) 



CC (tight) EC (tight) 



fiducial cuts 



shower shape 
isolation 
track match 



\A<j)\ > 0.02 
l^cail < 108 cm 

< 100 
/iso < 0.15 



IA0I > 0.02 — 

|2cai| < 108 cm 15 < < 25 

2trk| < 80 cm — 

< 100 < 100 

/iso < 0.15 /i,o < 0.15 

CTtrk < 5 (Jtrk < 10 



TABLE II. Number of W and Z candidate events. 

channel Z ee W ev 

fiducial region of electrons CC/CC CC/CC CC/EC CC 
e quality (t=tight, l=loose) t/1 t/t t/t t 
pass Main Ring Veto 537 1225 1116 28323 

fail Main Ring Veto 107 310 268 — 



FIG. 76. The difference between the predicted and the ac- 
tual 2-positions of the track center of gravity. 



TABLE III. Parton luminosity slope /3 and fraction of 
sea-sea interactions /ss in the W and Z production model. 
The j3 value is given for W ^ eu decays with the electron in 
CC and for Z ee decays with both electrons in CC. 







Z production 


W production 






13 (GeV-^) 


P (GeV-i) 


MRSA' [|5 




3.6 X 10"^ 


8.6 X 10"^ 0.207 


CTEQ3M 


42| 


3.3 X 10"^ 


8.7 X 10"^ 0.203 


CTEQ2M 


43) 




8.8 X 10"^ 0.203 


MRSD-' 1 




3.8 X 10"^ 


9.6 X 10"^ 0.201 
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TABLE IV. A'U|| for tlic W and Z event samples. 



Event Sample 


X;Six5/coshr,(e) (MeV) 


A'U|| (MeV) 


W 


95.8±0.4 


479±2±6 


Z 


93.6±1.3 


468±7±6 



TABLE V. Fitted values of gn for different parton distri- 
butitju fuiietioiis. L'licertaiiities are statistical oidv. 



prjee) < 15 GeV 



MRSA' 
MRSD-' 

CTEQ3M 
CTEQ2M 



0.59±0.10 GeV^ 
0.61±0.10 GeV^ 
0.54±0.10 GeV^ 
0.61±0.10 GeV^ 



A<f>(ee) 



0.64±0.14 GeV^ 
0.70±0.15 GeV^ 
0.57±0.13 GeV^ 
0.67±0.14 GeV^ 



TABLE VL The results of the Monte Carlo ensemble tests 
fitting the Mw mass for 105 samples of 28,323 events. 





mean 


rms 


correlation matrix 




(GeV) 


(GeV) 


rriT 


PT{e) 


Pt{v) 


rriT 


80.404 


0.067 


1 


0.669 


0.630 


Pr(e) 


80.415 


0.091 


0.669 


1 


0.180 




8().:!8I) 


0.105 


0.():!() 


0.180 


1 



TABLE IX. Uncertainties in the W mass measurement due 
to W production and decay model. 





rriT fit 


PT{e) fit 


Pt{v) fit 


Pt{W) spectrum 


10 MeV 


50 MeV 


25 MeV 


parton distribution functions 


20 MeV 


50 MeV 


30 MeV 


parton luminosity (3 


10 MeV 


10 MeV 


10 MeV 


radiative decays 


15 MeV 


15 MeV 


15 MeV 


W width 


10 MeV 


10 MeV 


10 MeV 


total 


30 MeV 


75 MeV 


45 MeV 



TABLE X. Variation of fitted W mass with choice of par- 
ton distribution function. 



rriT fit 



PT{e) fit 



Pxiv) fit 



MRSA' 
MRSD-' 
CTEQ3M 
CTEQ2M 





20 MeV 
5 MeV 
-21 MeV 





19 MeV 
48 MeV 
-17 MeV 





20 MeV 
22 MeV 
-30 MeV 



TABLE XL Changes in fitted W and Z masses if radiative 
effects are varied. 



variation 



rriT fit pT{e) fit pt(i^) fit m{ee) fit 



no radiative effects 50 MeV 43 MeV 
vary Ro byiO.l 3 MeV 4 MeV 



30 MeV 143 MeV 
MeV 19 MeV 



TABLE VII. The confidence levels from Kolmogorov- 
Smirnov tests comparing collider data to Monte Carlo pre- 
dictions for Mvi^ =80.44 GeV. 







rriT 


PT{e) 


Pt(i^) 


interval 


60- 


-90 GeV 


30-50 GeV 


30-50 GeV 


UT < 15 GeV 




0.25 


0.81 


0.20 


M|| < 




0.19 


0.78 


0.25 


W|| > 




0.61 


0.80 


0.48 


UT < 30 GeV 




0.55 


0.99 


0.58 


interval 


50- 


100 GeV 


25-55 GeV 


25-55 GeV 


UT < 15 GeV 




0.84 


0.83 


0.69 


W|| < 




0.77 


0.67 


0.62 


W|| > 




0.60 


0.66 


0.73 


UT < 30 GeV 




0.92 


0.80 


0.28 



TABLE XII. Uncertainties in the W mass measurement 
due to detector model parameters. 





mr fit 


prie) fit 


Pt{ 


ly) fit 


calorimeter linearity 


20 


MeV 


20 


MeV 


20 


MeV 


calorimeter uniformity 


10 


MeV 


10 


MeV 


10 


MeV 


electron resolution 


25 


MeV 


15 


MeV 


30 


MeV 


electron angle calibration 


30 


MeV 


30 


MeV 


30 


MeV 


electron removal 


15 


MeV 


15 


MeV 


20 


MeV 


selection bias 


5 


MeV 


10 


MeV 


20 


MeV 


recoil resolution 


25 


MeV 


10 


MeV 


90 


MeV 


recoil response 


20 


MeV 


15 


MeV 


45 


MeV 


total 


60 


MeV 


50 


MeV 


115 


MeV 



TABLE VIII. Uncertainties in the W mass measurement 
due to finite sample sizes. 





niT fit 


PT(e) fit 


priv) fit 


W sample 


70 MeV 


85 MeV 


105 MeV 


Z sample 


65 MeV 


65 MeV 


65 MeV 



total 95 MeV 105 MeV 125 MeV 



TABLE XIII. Uncertainties in the W mass measurement 
due to backgrounds. 

rriT fit pT{e) fit Pt{i') fit 

hadrons 10 MeV 15 MeV 20 MeV 

Z ^ee 5 MeV 10 MeV 5 MeV 

W ^ TV negligible 

cosmic rays negligible 

total 10 MeV 20 MeV 20 MeV 
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TABLE XIV. Summary of results from the 1992-1993 and 
1994-1995 data sets with the common and uncorrelated er- 
rors. 





1 QQ9 




1 QQ/1 


1 QQPi 


common 


Mw from rriT fit 


80.35 


GeV 


80.44 


GeV 






W statistics 


140 


MeV 


70 


MeV 






Z statistics 


160 


MeV 


65 


MeV 






calorimeter linearity 










20 


MeV 


calorimeter uniformity 










10 


MeV 


electron resolution 


70 


MeV 


20 


MeV 






electron angle calibration 










30 


MeV 


recoil resolution 


90 


MeV 


25 


MeV 






recoil response 


50 


MeV 


20 


MeV 






electron removal 


35 


MeV 


15 


MeV 






selection bias 


30 


MeV 


5 


MeV 






backgrounds 


35 


MeV 


10 


MeV 






W production/decay 










30 


MeV 


total uncertainty 


255 


MeV 


105 


MeV 


50 


MeV 



TABLE XV. Previously published measurements of the W 
boson mass. 



measurement 



Mw (GeV) 



reference 



CDF 90 
UA2 92 
CDF 95 
D0 96 
OPAL 96 
DELPHI 97 
L3 97 
ALEPH 97 



79.91±0.39 
80.36±0.37 
80.41±0.18 
80.35±0.27 

80.401^42 
80.40±0.45 
80.80lJ^:g 
80.14±0.35 
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